Gen. Rel. Grav., Dark Energy Special Issue manuscript No. 

(will be inserted by the editor) 



(N 



o 



X 



Received: 20 July 2007 / Accepted: 24 July 2007 



Thomas Buchert 

Dark Energy from structure: a status report 

S" 
o 
o 

(N 
o 

Q 

m . 

Abstract The effective evolution of an inhomogeneous universe model in any 
O I theory of gravitation may be described in terms of spatially averaged variables. 

<sJ '< In Einstein's theory, restricting attention to scalar variables, this evolution can be 

5^ ' modeled by solutions of a set of Friedmann equations for an effective volume scale 

bX). factor, with matter and backreaction source terms. The latter can be represented 

by an effective scalar field ('morphon field') modeling Dark Energy. 
', The present work provides an overview over the Dark Energy debate in con- 

^ ■ nection with the impact of inhomogeneities, and formulates strategies for a com- 

C<1 I prehensive quantitative evaluation of backreaction effects both in theoretical and 

■ observational cosmology. We recall the basic steps of a description of backreaction 
effects in relativistic cosmology that lead to refurnishing the standard cosmologi- 
cal equations, but also lay down a number of challenges and unresolved issues in 

' connection with their observational interpretation. 

I The present status of this subject is intermediate: we have a good qualitative 

understanding of backreaction effects pointing to a global instability of the stan- 
dard model of cosmology; exact solutions and perturbative results modeling this 
instability lie in the right sector to explain Dark Energy from inhomogeneities. It 
is fair to say that, even if backreaction effects turn out to be less important than 
anticipated by some researchers, the concordance high-precision cosmology, the 

■ architecture of current N-body simulations, as well as standard perturbative ap- 
proaches may all fall short in correctly describing the Late Universe. 
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1 General thoughts: 
— the standard model, the averaging problem and key insights 

1.1 Views on and beyond the standard model of cosmology 

The standard model of cosmology does not, like the standard model of particle 
physics, enjoy appreciable generality; it is based on the simplest conceivable class 
of (homogeneous-isotropic) solutions of Einstein's laws of gravitation. It is clear 
that the inhomogeneous properties of the Universe cannot be described by such a 
strong idealization. The key issue is whether they can be described so on average, 
and this is the subject of considerable debate and controversy in the recent litera- 
ture. If the standard model indeed describes the averaged model, we have to show 
that backreaction effects, being the main subject of this report, are negligible. We 
are striving to discuss most of the related aspects of this debate. 

1.1.1 Dark Energy and Dark Matter 

In the standard model of cosmology one has to conjecture the existence of two 
constituents, if observational constraints are met, that both have yet unknown ori- 
gin: first, a dominant repulsive component is thought to exist that can be mod- 
eled either by a positive cosmological constant or a scalar field, e.g. a so-called 
quintessence field. Besides this Dark Energy, there is, secondly, a non-baryonic 
component that should considerably exceed the contribution by luminous and dark 
baryons and massive neutrinos. This Dark Matter is thought to be provided by 
exotic forms of matter, not yet detected in (non-gravitational) experiments. Ac- 
cording to the concordance model fTTS\, (Fl, flSO |, the former converges to about 
3/4 and the latter to about 1 /4 of the total source of Friedmann's equations, up to 
a few percent that have to be attributed to baryonic matter and neutrinos (in the 
matter-dominated era). There are, however, other voices [19j . ifTSl . 

Contemporary research to uncover this enigma pursues essentially two direc- 
tions: one focusses on generalizations of the geometry of spacetime mostly re- 
stricting attention to modifications of the underlying theory of gravitation, the 
other invokes new sources in the energy momentum tensor and so implies a chal- 
lenge for particle physics. As for the former, a Dark Energy component may pos- 
sibly derive either from higher-order Ricci curvature Lagrangians II53I (as well as 
Capozziello and Francaviglia, this volume), [67], or string-motivated low-energy 
effective actions |20|. It is doubtful whether a fundamental scalar field exists in 
nature, at least one that can be viewed as a natural candidate for the relevant ef- 
fects needed to explain Dark Energy. This latter remark is supported by the well- 
known violation of energy conditions of a quintessence field that is able to pro- 
duce late-time volume acceleration of the Universe. Rather, a scalar field would 
likely be an effective one, either stemming from higher-order gravity terms, or 
effective terms as remnants from higher dimensions that are compactified or even 
non-compactified as in brane world cosmologies 111271 (see also Koyama, this 
volume). As we shall learn below, already classical general relativity allows to 
identify effective geometrical terms, simply resulting from inhomogeneities, with 
an effective scalar field component, the morphon field |48|, a good example of 
William of Ockham's razor. In this picture Dark Energy emerges as an excess of 
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kinetic over potential energies of a scalar field in an 'out-of-equilibrium' state, 
and it allows attributing Dark Energy to the classical vacuum. If we restrict our 
attention to cosmology and the fitting of extra terms from various different mod- 
ified gravitational theories to observational data, then those extra terms may also 
be mapped into morphon fields with different but unambiguously defined physical 
consequences. A review of the status and properties of currently discussed mod- 
els can be found in ll65l . see also II143II (as well as Padmanabhan, this volume), 
II189I , 0167 1. We shall not directly address the Dark Matter problem in this report, 
but also this problem might be related to an explanation of Dark Energy; we shall 
discuss such possible relations. 

Thus, the intriguing question is whether an explanation of these dark com- 
ponents is (i) the task of particle physicists, or (ii) an expression of the need to 
modify the laws of gravitation, or (iii) whether the cosmological model is built on 
oversimplified priors. We are going to study this last possibility. 

1.1.2 The longstanding averaging problem 

Does an inhomogeneous model of the Universe evolve on average like a homoge- 
neous solution of Einstein's or Newton's laws of gravitation? This question is not 
new, at least among relativists who think that the answer is certainly, in general, 
no, not only in view of the nonlinearity of the theories mentioned II70I . The prob- 
lem was and still is the notion of averaging whose specification and unambiguous 
definition turned out to be an endeavor of high magnitude, mainly because it is not 
straightforward to give a unique meaning to the averaging of tensors, e.g., a given 
metric of spacetime. This problem seems to lie in the backyard of relativists who, 
from time to time, add another effort towards a solution of this technical issue. On 
the other hand, the community of cosmologists should locate exactly this research 
topic at the basis of their evolutionary models of the Universe. 

Although there have been numerous exceptions to this ubiquitous ignorance 
of the averaging problem in cosmology, e.g. II177II . and many efforts after George 
Ellis ITOJ has brought the subject into the fore, ESI, El, 113, lfT05l . i56|, I86J, 
lfT96l , 1761 . lfT68]| . IITI . ifTPfT. still, the cosmologist's thinking rests on the hege- 
mony of the standard model despite the drastic changes of our picture of struc- 
tures in the Universe on large scales. This standard model, up to the present state 
of knowledge, is used as a prior to interpret a wide variety of orthogonal obser- 
vations, and it is therefore hard to beat due to this intentionally established status. 
Therefore, most investigations in cosmology are still based on the vocabulary of 
the standard model, aiming to constrain its global cosmological parameters, often 
on the basis of observations of structure in the regional Universe that is very dif- 
ferent from homogeneous and isotropic. As a consequence, also structure on large 
scales is described in terms of (quasi-Newtonian) perturbations of this standard 
model, a construction that again makes only sense, if the standard model correctly 
describes the average distributions of matter and geometry. Promisingly, the con- 
jecture that the standard model agrees with the averaged model has recently been 
recognized as such and challenged by a wider community thanks to the Dark En- 
ergy debate. 



This is certainly an incomplete list - more references may be found in these papers and, 
e.g., in l72i. 
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1.1.3 Uncharted territory beyond the standard model 

The concordance model is encircled by a large set of observational data that 
are, however, orthogonal only within the predefined solution space of a FLRW 
(Friedmann-Lemaitre-Robertson- Walker) cosmology. This solution space has di- 
mension two for Friedmann's expansion law derives from the Hamiltonian con- 
straint of general relativity (see Eq. ( 118b below), restricted to (about every point) 
locally isotropic and hence (by Schur's Lemma) homogeneous distributions of 
matter and curvature, 

n,n+nk + ^A = 1 , (1) 

where the standard cosmological parameters are global and iconized by the cosmic 
triangle lllil . 

pH (f ) is the homogeneous matter density, H{t) ■.= a/a Hubble's function with the 
scale factor a(f), k a positive, negative or vanishing constant related to the three 
elementary constant-curvature geometries, and A is the cosmological constant, 
nowadays - if positive - employed as the simplest model of Dark Energy II151I . 

We shall learn below that an extended solution space of an averaged inhomo- 
geneous universe model is three-dimensional, when we include inhomogeneities 
of matter and geometry. Hence, such more realistic models seem to enjoy more pa- 
rameter freedom, but it should be emphasized that these (effective) 'parameters' 
are defined in terms of volume averages of dynamically interacting physical vari- 
ables. For a given inhomogeneous model, the additional parametrization appears 
in the initial conditions for the inhomogeneities that are absent in the standard 
model of cosmology. 

How can we be sure that fitting an idealized model, that ignores inhomo- 
geneities, to observational data is not 'epicyclic', especially if the model enters 
as a prior into the process of interpreting the data? Confronting observers with the 
wider class of averaged cosmologies allows them to draw their data points within 
a cube of possible solutions and to differentiate the relevant observational scales 
reflected by these data; if we 'force' them to draw the data points into the plane 
of the FLRW solutions on every scale, then they conclude that there are 'dark' 
components. Thus, we have to exclude that they may have missed something in 
the projection and we have to clarify whether the ignorance of scale-dependence 
of observables in the standard model does not mislead their interpretation. Both 
issues are equally important to judge the viability of the standard model in ob- 
servational cosmology: the first is the question of how backreaction quantitatively 
affects the standard cosmological parameters, and the second is the comparison of 
data taken on small scales (e.g. on cluster scales) and data taken on large scales 
(e.g. CMB; high-redshift supernovae). Both additional 'degrees of freedom' in 
interpreting observational data are interlocked in the sense that backreaction ef- 
fects may alter the evolution history of cosmological parameters. A comparison 
of data taken on different spatial scales has therefore also to be subjected to a crit- 
ical assessment of data that are taken at different times of the cosmic history: with 
backreaction at work, the simple time-scaling of parameters in a FLRW cosmol- 
ogy is also lost. 
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The plan of this report is the following. We shall first provide a list of argu- 
ments that justify existence of backreaction effects. Then, we move on to construct 
realistic universe models and discuss the governing equations in Section 2. A qual- 
itative understanding of the backreaction mechanism relevant to the question of 
Dark Energy is developed in Section 3, and thereafter we propose and discuss 
strategies for a quantitative evaluation of backreaction effects in Section 4. Before 
we now enter the physics of backreaction that is easy to understand, we have to 
probe some more critical territory in the following subsection. 

1.2 Averaging strategies: different 'directions' of backreaction 

The notion of averaging in cosmology is tied to space-plus-time thinking. Despite 
the success of general covariance in the four-dimensional formulation of classical 
relativity, the cosmologist's way of conceiving the Universe is evolutionary. This 
breaking of general covariance is in itself an obstacle to appreciating the proper 
status of cosmological equations. The standard model of cosmology is employed 
with the implicit understanding that there is a global spatial frame of reference 
that, if mapped to the highly isotropic Cosmic Microwave Background, is elevated 
to a physical frame rather than a particular choice of a mathematical slicing of 
spacetime. Restricting attention to an irrotational cosmic continuum of dust (that 
we shall retain throughout the main text), the best we can say is that all elements 
of the cosmic continuum defined by the homogeneous distribution of matter are in 
free fall within that spacetime, and therefore are preferred relative to accelerating 
observers with respect to this frame of reference. Those preferred observers are 
called fundamental. Exploiting the diffeomorphism degrees of freedom we can 
write the FLRW cosmology in contrived ways, so that nobody would realize it 
as such. This point is raised as a criticism of an averaging framework |99|, as 
if this problem were not there in the standard model of cosmology. Again, the 
'natural' choice for the matter model 'irrotational dust' is a collection of freely- 
falling continuum elements, now for an inhomogeneous continuum. For such a 
generalized collection of fundamental observers, the 4-metric form readfl 

4g = -dt^ +'g-'g= gab dX" ® dX" , (3) 

where latin indices run through 1 • • • 3 and X" are local (Gaussian normal) coor- 
dinates. Evolving the first fundamental form of the spatial hypersurfaces along 
d/dt =: dt defines their second fundamental form 

= dX" ® dX" ■ Kab:= -^drgab , (4) 

with the extrinsic curvature components Kaj,. Such a comoving (synchronous) slic- 
ing of spacetime may be considered 'natural', but it may also be questioned. How- 
ever, to dismiss its physical relevance due to the fact that shell-crossing singular- 
ities arise is shortsighted. It is a problem of the matter model in the first place. A 

' For notations the reader may consult the Appendix; generally, we work with spatial variables 
in the hypersurfaces of constant coordinate time t (that is equal to proper time for an iiTotational 
dust continuum), and we explicitly indicate with a prefix when we talk about four-dimensional 
variables in cases where this is not obvious. 
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comoving (Lagrangian) frame helps to access nonlinear stages of structure evo- 
lution, as is well-exemplified in Newtonian models of structure formation, where 
the problem of choosing a proper slicing is absent. Those nonlinear stages in- 
evitably include the development of singularities, provided we do not improve on 
the matter model to include effects that counteract gravitation (like velocity dis- 
persion) in order to regularize such singularities |43 1 . If a chosen slicing appears 
to be better suited, because it does not run into singularities, then one should rather 
ask the question whether the evolution of variables is restricted to a singularity- 
free regime just because inhomogeneities are not allowed to enter nonlinear stages 
of structure evolution. An example for this is perturbation theory formulated e.g. 
in longitudinal gauge, where the variables are 'gauge-fixed' to a (up to a given 
time-dependent scale factor) non-evolving background. 

However, the problem of choosing an appropriate slicing of spacetime is not 
off the table. There exist strategies to consolidate the notion of an effective spa- 
tial slicing that would minimize frame fluctuations being attributed to the diffeo- 
morphism degrees of freedom in an inhomogeneous model. Such, more involved, 
strategies relate to the intrinsic direction of backreaction that we put into perspec- 
tive below. 

1.2.1 Extrinsic (kinematical) and intrinsic backreaction 

Having chosen foliation of spacetime implies that we can speak of two 'direc- 
tions' : one being extrinsic in the direction of the extrinsic curvature Kah of the em- 
bedding of the hypersurface into spacetime (e.g. parametrized by time), the other 
being intrinsic in the direction of the Ricci tensor Rgi, of the three-dimensional 
spatial hypersurfaces parametrized by a scaling parameter (let it be the geodesic 
radius of a randomly placed geodesic ball). Consequently, we may speak of two 
'directions' of backreaction: inhomogeneities in extrinsic curvature and in intrin- 
sic curvature. The former is of kinematical nature, since we may interpret the ex- 
trinsic curvature actively through the expansion tensor Qat '■= —Kah, and introduce 
a split into its kinematical parts: 0„ft = 1 /3g„h0 + Oat, with the rate of expansion 
= 0^,, the shear tensor Gat, and the rate of shear := X/lOabO"^; note that 
vorticity and acceleration are absent for dust in the present flow-orthogonal foli- 
ation. The latter addresses the so-called ^fring problem [701, [761, |fl48|, i.e. the 
question whether we could find an effective constant-curvature geometry that best 
replaces the inhomogeneous hypersurface at a given time. An answer to this ques- 
tion has to deal with the problem of 'averaging' the tensorial (spatial) geometry for 
which several different strategies are conceivable. Some of those strategies do not 
distinguish between extrinsic and intrinsic averaging (e.g. |196|, [62], [63l, and 
other references in |72|). A comparison of such a more 'synthetic' approach with 
a pure kinematical averaging that leaves the physical properties of a spatial hyper- 
surface untouched has been provided 1 149] and helps to also formally understand 
the differences between both viewpoints. 

One method has recently obtained a strong position in the context of Perel- 
man's work (e.g. II154II . 111558 ) on the Ricci-Hamilton flow related to the recent 
proof of Poincare's conjecture, and implied progress on Thurston's geometrization 
program |5 | to cut a Riemannian manifold into 'nice pieces' of eight elementary 
geometries. This method we briefly sketch now. 
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1.2.2 Renormalization of average characteristics: smoothing the geometry 

Employing the Ricci-Hamilton flow f911, \92\, \55\, an 'averaging' of geometry 
can be put into practice by a rescaling of the spatial metric tensor, much in the 
spirit of a renormalization flow 11561 . A general scaling flow is described by Pe- 
tersen's equations 1 156] that we may implement through a 2+1 setting by evolving 
the boundary of a geodesic ball in a three-dimensional cosmological hypersurface 
in radial directions, thus exploring the Riemannian manifold passively. Upon lin- 
earizing the general scaling flow, e.g. in normal geodesic coordinates, we obtain a 
scaling equation for the metric along radial directions; up to tangential geometrical 
terms on the boundary we obtain II39I , 

§^8ab{r) - jzgab{r) = -2Rab{rQ)[r ~ ro] , (5) 
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i.e. the metric scales in the direction of its Ricci tensor much in the same way 
as it is deformed in the direction of the extrinsic curvature by the Einstein flow. 
If we now implement the active (geometrically Lagrangian) point of view of de- 
forming the metric by the same flow along a Lagrangian vector field d / dro while 
holding the geodesic radius tq fixed, we are able to smooth the metric in a con- 
trolled way. Depending on our choice of normalization of the flow, we may pre- 
serve the mass content inside the geodesic ball while smoothing the metric. Such 
a mass-preserving Ricci flow transforms kinematical averages on given hypersur- 
faces from their values in the inhomogeneous geometry (the actual space section) 
to their values on a constant-curvature geometry (the fitting template for the space 
section): they are renormalized resulting in additional backreaction effects due to 
the difference of the two volumes (the Riemannian volume of the actual space sec- 
tion and the constant-curvature volume) - the volume effect, and also curvature 
backreaction terms that involve averaged invariants of the Ricci tensor. For details 
and references see (391 and for small overviews |40| and [41]. In such a setting 
the role of lapse and shift functions (i.e. the choice of slicing, cf. Appendix) can 
also be controlled by employing the recent results of Perelman li54l . 

We now come to some crucial points of understanding the physics behind 
backreaction. In order not to think of any exotic mechanism, the historical use 
of the notion 'model with backreaction' should simply be replaced by 'more real- 
istic model' . 



1.3 The origin of kinematical backreaction and the physics behind it 

Let us now concentrate on the question, why there must be backreaction at work, 
restricting attention to kinematical backreaction as defined above. In doing so, 
we do not actively modify the physics, i.e. the metrical properties of spatial sec- 
tions; we merely look at general integral properties of the inhomogeneous spatial 
distributions of matter and geometry on a given scale. After we have understood 
the reasons behind backreaction effects in general terms, i.e. without resorting to 
restrictions of spatial symmetry or approximations of evolution models, the very 
question of their relevance is better defined. 
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1.3.1 An incomplete message to particle physicists 

Employing Einstein's general theory of relativity to describe the evolution of the 
Universe, we base our universe model on a relation between geometry and mat- 
ter sources. A maximal reduction of this theoretical fundament is to consider the 
simplest conceivable geometry. Without putting in doubt that it might be an over- 
simplification to assume a (about every point) locally isotropic (and hence homo- 
geneous) geometry, standard cosmology conjectures the existence of sources that 
would generate this simple geometry. As already remarked, the majority of these 
sources have yet unknown physical origin. Obviously, particle physicists take the 
demand for missing fundamental fields literally. But, as was emphasized above, 
the standard model has physical sense only, if a homogeneous-isotropic solution 
of Einstein's equations also describes the inhomogeneous Universe effectively, i.e. 
on average. This is not obvious. The very fact that the distributions of matter and 
geometry are inhomogeneous gives rise to backreaction terms; we shall restrict 
them to those additional terms that influence the kinematics of the homogeneous- 
isotropic solutions. These terms can be viewed to arise on the geometrical side of 
Einstein's equations, but they may as well be put on the side of the sources. 

We start with a basic kinematical observation that lies at the heart of the back- 
reaction problem. 

1.3.2 A key to the averaging problem: non-commutativity 

Let us define spatial averaging of a scalar field f on a compact domain ^ with 
volume := \ ^\ through its Riemannian volume average 

{^{X\t))^:=^I^W{X\t)Jd'X ; /:=^detfe,). (6) 

The key property of inhomogeneity of the field W is revealed by the commutation 
rule mi, W 

d,{W)r,^{d,W)^ = {0W)^~{0)^{W)^ , (7) 

where := m^^ denotes the trace of the fluid's expansion tensor, u^ its 4-velocity, 
and dtJ = 0J the evolution of the root of the 3-metric determinant 7; the spatial 
average of describes the rate of volume change of a collection of fluid elements 
along d I dt, 

{0)^ = ^=:3H^, (8) 

where we have introduced a volume Hubble rate that reduces to Hubble's 
function in the homogeneous case. Commutativity reflects the conjecture implied 
by the standard model: a realistically evolved inhomogeneous field will feature 

^ This is a strong assumption on smaller spatial scales in the case of the matter model 'irrota- 
tional dust': as soon as singularities in the flow develop, the boundary of the domain then also 
experiences singularities, i.e. a breaking of the boundary due to a splitting of the domain or due 
to a merging of domains. These latter processes that alter the domain's topology may also occur 
in a smooth way, if the flow is regularized through generalizations of the matter model. 
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the same average characteristics as those predicted by the evolution of the (homo- 
geneous) average quantity; in other words, the right-hand-side of (|7]) is assumed 
to vanish. This rule also shows that backreaction terms deal with the sources of 
non-commutativity that are in general non-zero for inhomogeneous fields. Note 
that this rule is purely kinematical, which shows that it is not necessarily the non- 
linearity of the field equations that is responsible for backreaction effects. 

1.3.3 Regional volume acceleration despite local deceleration 

Based on a first application of the above rule, we shall emphasize that there is 
not necessarily anti-gravity at work, e.g. in the 'redcapped' version of a positive 
cosmological constant, in order to have sources that counteract gravity. Raychaud- 
huri's equation, if physically essential terms like vorticity, velocity dispersion, or 
pressure are retained, provides terms needed to oppose gravity, e.g., to support 
spiral galaxies (vorticity), elliptical galaxies (velocity dispersion), and other sta- 
bilization mechanisms involving pressure (think of the hierarchy of stable states 
of stars until they collapse into a Black Hole). Admittedly, those terms are effec- 
tively 'small-scale-players'. Now, let us consider Raychaudhuri's equation (see 
([2TI below), restricted to irrotational dusQ, 

dfS =A-AnGp + 2ll-l^ , (9) 

with the principal scalar invariants of @„i„ 111 : = 2/30^ — 2a^ and I : = 0. Then, 
unless there is a positive cosmological constant, there is no term that could counter- 
balance gravitational attraction and, at every point, (9,0 < 0. Applying the com- 
mutation rule (ITI for f = 0, we find that the averaged variables obey the same 
equation as above despite non-commutativitj0: 

d,{e)^=A-AnG{p)^ + 2{ll)^-{l)l . (10) 

This result can be understood on the grounds that shrinking the domain ^ to a 
point should produce the corresponding local equation. Now, notwithstanding, the 
above equation contains a positive term that acts against gravity. This can be easily 
seen by rewriting the averaged principal invariants: we obtairO 

2(11)^- (1)2, = ^((0-(0)^)2)^-2<(cT-(a>^)2)^-i(0>|-2(cT)2^ , 

(11) 

which, compared with the corresponding local expression, 

211-1^ = -^02-2a2 , (12) 

^ We assume that the influence of a strong vorticity evolution (that is known to happen on 
small scales in the nonlinear regime of structure formation) is not relevant on scales larger than 
the scale of, say, superclusters of galaxies. According to the sign of its appearence in Raychaud- 
huri's equation, vorticity counteracts gravitation and its effect will be relevant, if averages are 
performed over domains on and below the scales of galaxy clusters. 

^ This is only true, if all terms appearing in Raychaudhuri's equation are written in terms of 
principal scalar invariants; it is actually a special non-linearity of this equation that cancels the 
corresponding non-commutativity term (see Corollary I in 1 32 1). 

^ We have formally inserted the averaged shear term, so that the last two terms correspond to 
the local ones. 
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gave rise to two additional, positive-definite fluctuation terms, where that for the 
averaged expansion variance enters with a positive sign. It may appear 'magic' 
that the time-derivative of a (on some spatial domain ^) averaged expansion may 
be positive despite the fact that the time-derivative of the expansion at all points 
in S> is negative. As the above explicit calculation shows, this property does not 
furnish an argument against the possibility of volume acceleration [991 , but simply 
is due to the fact that an average correlates the local contributions, and it is this 
correlation (or fluctuation) that adds 'kinematical pressure' . The interesting point 
is that these additional terms are 'large-scale players', as we shall make more 
precise belowQ. 

What we can learn from this simple exercise is that any local argument, e.g. 
on the smallness of some perturbation amplitude at a given point, is not enough 
to exclude regional ('global') physical effects that arise from averaging inhomo- 
geneities; even if deviations from the average are small, as measured for example 
today, the evolution of the average may be different from the evolution of a 'back- 
ground solution' in perturbation theory. As we shall discuss more in detail in the 
course of this report, such correlation effects must not be subdominant compared 
to the magnitude of the local fields, since they are related to the spatial variation 
of the local fields and, having said 'spatial', it could (and it will) imply a coupling 
to the geometry as a dynamical variable in Einstein gravitation. This latter remark 
will turn out very useful in understanding the potential relevance of backreaction 
effects in relativistic cosmology. 

1.3.4 The production of information in the Universe 

The above considerations on effective expansion properties can be essentially 
traced back to 'non-commutativity' of averaging and time-evolution, lying at the 
root of backreaction. (Note that additional 'spatial' backreaction terms that have 
been discussed in Subsect. [l.2.2l are also the result of a 'non-commutativity', this 
time between averaging and spatial rescaling - see also II72II .') The same reasoning 
underlies the following entropy argument. Applying the commutation rule (|7]l to 
the density field, W = p, 



we derive, as a source of non-commutativity, the (for positive-definite density) 
positive-definite Lyapunov functional (known as Kullback-Leibler functional in 
information theory; |96| and references therein): 



This measure vanishes for Friedmannian cosmologies ('zero structure'). It attains 
some positive time-dependent value otherwise. The source in ( 113b shows that rel- 
ative entropy production and volume evolution are competing: commutativity can 

^ The physical and observational consequences of the expansion fluctuation term have been 
thoroughly explained and illustrated by a toy model in the review paper 116111 . 




(13) 




(14) 
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be reached, if the volume expansion is faster than the production of information 
contained within the same volume. 

In [|96i the following conjecture was advanced: 

The relative information entropy of a dust matter model S{p\ \ (p)^} is, for 
sufficiently large times, globally (i.e. averaged over the whole manifold L that is 
assumed simply-connected and without boundary) an increasing function of time. 

This conjecture already holds for linearized scalar perturbations at a Friedmannian 
background (the growing-mode solution of the linear theory of gravitational insta- 
bility implies d,S>0 and S is, in general, time-convex, i.e. d^S > 0). Generally, 
information entropy is produced, i.e. dfS > with 

^'^^P^l^P^^^ = {5p0)r, = - (p50)^ = - {5p50)^ , (15) 

(and with the deviations of the local fields from their average values, e.g. 5p : = 
p — (p) if the domain S> contains more expanding underdense and contracting 
overdense regions than the opposite states contracting underdense and expanding 
overdense regions. The former states are clearly favoured in the course of evolu- 
tion, as can be seen in simulations of large-scale structure. 

There are essentially three lessons relevant to the origin of backreaction that 
can be learned here. First, structure formation (or 'information' contained in struc- 
tures) installs a positive-definite functional as a potential to increase the devi- 
ations from commutativity; it can therefore not be statistically 'averaged away' 
(the same remark applies to the averaged variance of the expansion rate discussed 
before). Second, gravitational instability acts in the form of a negative feedback 
that enhances structure (or 'information'), i.e. it favours contracting clusters and 
expanding voids. This tendency is opposite to the thermodynamical interpreta- 
tion within a closed system where such a relative entropy would decrease and the 
system would tend to thermodynamical equilibrium. This is a result of the long- 
ranged nature of gravitation: the system contained within & must be treated as an 
open system. Third, backreaction is a genuinely non-equilibrium phenomenon, 
thus, opening this subject also to the language of non-equilibrium thermodynam- 
ics Us?], [176], [200], general questions of gravitational entropy II152I , II153I . 
II24I . II96I . II135I . and observational measures using distances to equilibrium II13I . 
'Near-equilibrium' can only be maintained (not established) by a simultaneous 
strong volume expansion of the system. Later we discuss an example of a cosmos 
that is 'out-of-equilibrium', i.e. settled in a state far from a Friedmannian model 
that, this latter, can be associated with the relative equilibrium state 5 = 0. 

In particular, we conclude that the standard model may be a good descrip- 
tion for the averaged variables only when information entropy production is over- 
compensated by volume expansion (measured in terms of a corresponding adi- 
mensional quantity). This latter property is realized by linear perturbations at a 
FLRW background. Thus, the question is whether this remains true in the non- 
linear regime, where information production is strongly promoted by structure 
formation and expected to be more efficient. 

Before we can go deeper into the problem of whether such backreaction terms, 
being well-motivated, are indeed relevant in a quantitative sense, we have to study 
the governing equations. 
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2 Constructing a realistic universe model: 
— refurnishing the cosmological equations 

In this section we recall a set of averaged Einstein equations together with alter- 
native forms of these equations which put us in the position to study backreaction 
terms as additional sources to the standard Friedmann equations. 

2.1 Einstein's equations recalled 

In order to make the presentation more self-contained, we recall the complete set 
of local Einstein equations, restricted to irrotational fluid motion with the sim- 
plest matter model 'dust' (i.e. vanishing pressure), as befor^l. In this case the 
flow is geodesic and space-like hypersurfaces can be constructed that are flow- 
orthogonal at every spacetime event in a 3 + 1 representation. 
We start with Einstein's equation^ 

'^Rfiv -^8liv'^R = 8^GpUfiUv -Ag^v , (16) 

with the 4-Ricci tensor ^R^v, its trace ^/?, the fluid's 4- velocity (u^u^ = - 1), 
the cosmological constant A , and the rest mass density p obeying the conservation 
law 

ipu^u").^ = . (17) 

In a flow-orthogonal coordinate system = {X'^j) (i.e., Gaussian or normal 
coordinates which are comoving with the fluid) we can write — f^{X^,t), and 
we have =/'' = ( 1 , 0, 0, 0) and = = (— 1 , 0, 0, 0) . These coordinates are 
defined such as to label geodesies in spacetime, i.e., u^u^.^ = 0. 

Defining the two fundamental forms as in Eqs. QHI, with the 3-metric coeffi- 
cients gij and the extrinsic curvature coefficients Kij := —h^-h'^jU^-y (projected into 
the hypersurfaces orthogonal to with the help ofh^y : = ^/j v + "/j "v)^ Einstein's 
equations ( fT6] l together with ( fTTt (contracted with Uy) then are equivalent to the 
following system of equations fJ], fTTE\, consisting of the energy or Hamiltonian 
constraint and the momentum or Codazzi constraints, 

^(r + K^- k'jK^^ =SnGp+A ; K)^ , . -Kij=0 , (18) 

and the evolution equations for the density and the two fundamental forms, 

drp = Kp ; drg^j = -2 g^kK) ; d^K' j = KK'j+R<j - {AnGp +A)5) . (19) 

* The corresponding equations with arbitrary lapse and shift functions for a perfect fluid 
energy-momentum-tensor are discussed in the Appendix, together with the averaged equations. 

' Greek indices run through 0...3, while latin indices run through 1...3; summation over re- 
peated indices is understood. A semicolon will denote covariant derivative with respect to the 
4-metric with signature (—,-!-, the units are such that c = 1; further below, a double ver- 

tical slash II denotes covariant derivative with respect to the 3— metric g,,, while a single vertical 
slash denotes partial derivative with respect to the local coordinates X'\ The overdot denotes 
partial time-derivative (at constant X') as before, here identical to the covariant time-derivative 
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R '■= R'j and K : = K\ denote the traces of the spatial Ricci tensor Rij and the 
extrinsic curvature Kij, respectively. Expressing the latter in terms of kinematical 
quantities, 

- Kij = 0,j = Gij + ^-0gij ■ ~K = , (20) 

with the expansion 0,,, the trace-free symmetric shear a,y , and the rate of expan- 
sion 0, we may write the above equations in the form 

^R + ^0^-a^ = 8nGp+A ; a^||. = ^0|- 

2 

d,p = 0p ■ d,gij = 2 giko) + -0gik5) ; 
d,0 + ^0^ + 2a^ + 4nGp-A = ; 

d,a) + 0a) = -I^R'j-\d)R^ , (21) 

where we have introduced the rate of shear := l/2a'^a'. (To derive the last 
two equations, Raychaudhuri 's equation fl631 and the equation for the trace-free 
parts, we have used the Hamiltonian constraint.) 



2.2 Averaged cosmological equations 

In order to find evolution equations for effective (i.e. spatially averaged) cosmo- 
logical variables, we may put the following simple idea into practice. We observe 
that Friedmann's differential equations ll82l . ll83ll capture the scalar parts of Ein- 
stein's equations ( 1211 ). while restricting them by the strong symmetry assump- 
tion of local isotropy. The resulting equations, Friedmann's expansion law (the 
energy or Hamiltonian constraint) and Friedmann's acceleration law (Raychaud- 
huri's equation), together with restmass conservation, 

\ 2 / 

d\ „ . 3k , „ . ^ I d 



3 - -StzGPh-A = ; 3-+4nGpH-A = ; P/f+3 - ]Ph=0 , 

(22) 

can be replaced by their spatially averaged, general counterparts (for the details 
the reader is referred to lf321[34l[3^|48l '): 

3('-^\ ~^nG{p),-A = ^^^1^ ; (23) 



a^ 



3^+4nG{p)^~A = ; (24) 



{ph + 3^{p)^ = . (25) 
a® 



14 



We have replaced the Friedmannian scale factor by the volume scale factor a^, 
depending on content, shape and position of the domain of averaging f^, defined 
via the domain's volume 'V^{t) = \^\, and the initial volume V^. = Vsi{ti) = |^,|: 

-<')-(^)"^ 

Using a scale factor instead of the volume should not be confused with 'isotropy' . 
The above equations are general for the evolution of a mass-preserving, com- 
pact domain containing an irrotational continuum of dust, i.e. they provide a 
background-free and non-perturbative description of inhomogeneous and aniso- 
tropic 

fieldS 

The new term appearing in these equations, the kinematical back- 
reaction, arises as a result of expansion and shear fluctuations: 



:= 2 (11)^ - - {1)1 = 3 ((0 - {0)^f)^ - 2(ct2>5, ; (27) 

I and II denote the principal scalar invariants of the extrinsic curvature, and the 
second equality follows by introducing the decomposition of the extrinsic curva- 
ture into the kinematical variables, as before. Also, it is not a surprise that the 
general averaged 3-Ricci curvature (R) ^ replaces the constant-curvature term in 
Friedmann's equations. Note also that the term encoding the fluctuations has 
the particular structure of vanishing at a Friedmannian background, a property that 
it shares with gauge-invariant variable^. 

In the Friedmannian case, Eqs. (122b . the acceleration law arises as the time- 
derivative of the expansion law, if the integrability condition of restmass conser- 
vation is respected, i.e. the homogeneous density 'x . In the general case, 
however, restmass conservation is not sufficient. In addition to the (built-in) gen- 
eral integral of Eq. ( 1251 ), 



( p(f/))g, _ 



{P)9 = '-^—^ = ; = M^, , (28) 



we also have to respect the following curvature-fluctuation-coupling: 

-^dJQ^a%) + -Ld,{{R)^al) = 0- (29) 
a^ ^ ^ 

This relation will be key to understand how backreaction can take the role of Dark 
Energy. 

One could, of course, introduce an isotropic or anisotropic reference background 1441 or, 
explicitly isolate an averaged shear from the above equations to study deviations from the kine- 
matics of Bianchi-type models, as was done with some interesting results in 1 12|. 

" In a quasi-Newtonian setting, where averages are taken on the Euclidean or constant- 
curvature background space, the variable Qjs is gauge-invariant to se cond- o rder in perturba- 
tion theory 1 1 16|, 1 123 1, since this variable vanishes at the background 11811 , 11821 ; for related 
thoughts see L149J . LI 46 J . 
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2.3 Alternative forms of the averaged equations 

We here provide three compact forms of the averaged equations introduced above, 
as well as some derived quantities. They will prove useful for our further discus- 
sion of the backreaction problem. 

2.3.1 Generalized expansion law 

The correspondence between Friedmann's expansion law (the first equation in 
(122b ') and the general expansion law ( 1231) can be made more explicit through formal 
integration of the integrability condition (129b : 



The (domain-dependent) integration constant ke^. relates the new terms to the 
'constant-curvature part'. We insert this latter integral back into the expansion 
law ( 123b and obtain: 



This equation is formally equivalent to its Newtonian counterpart ll44ll . It shows 
that, by eliminating the averaged scalar curvature, the whole history of the aver- 
aged kinematical fluctuations acts as a source of a generalized expansion law that 
features the 'Friedmannian part' on the left-hand-side of ( 1311) . 

2.3.2 Effective Friedmannian framework 

We may also recast the general equations ( 1231 l24l l25l |29] | by appealing to the 
Friedmannian framework. This amounts to re-interpret geometrical terms, that 
arise through averaging, as effective sources within a Friedmannian setting. 

In the present case the averaged equations may be written as standard zero- 
curvature Friedmann equations for an effective perfect fluid energy momentum 
tensor with new effective sources 1 34 1 : 




(30) 




(31) 



Pel = (P) 



1 



1 




\6kG 
1 



\6nG 
1 




l6nG 



Q9 + 



487rG 




(32) 




3^+47rG(pl + 3pfff)-A 



; 




(33) 
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Eqs. ([33]) correspond to the equations (|23] |, (|24] |. ( |25]) and ( |29]) , respectively. 

We notice that 2;^, if interpreted as a source, introduces a component with 
'stiff equation of state', Pq = Pg, suggesting a correspondence with a free scalar 
field (discussed in the next subsection), while the averaged scalar curvature intro- 
duces a component with 'curvature equation of state' = — l/3p^. Although 
we are dealing with dust matter, we appreciate a 'geometrical pressure' in the 
effective energy-momentum tensor. 

There is, of course, some ambiguity in defining the effective sources. We recall 
||36J that, firstly, it may sometimes be useful to incorporate A into the effective 
sources by defining p^j^ ■= P^ + A/StzG and pffi/^ ■= pff^ — A/SnG. Secondly, 
we might add the 'constant-curvature term' 3k^./ajg to the expansion law in ( [33] l; 
if we wish to do so, then the effective sources can be represented solely through 
the kinematical backreaction term and its time-integral. For this we have to 
exploit the 'Newtonian form', Eq. ( 131b , and would have to define the effective 
sources as follows: 

(34) 

The integrated form of the integrability condition, Eq. ( l30l) , then allows to express 
again through the averaged scalar curvature, Xq, = (ike^.jc^^ — Qgj — (/?)^, and 
we obtain the sources corresponding to ( 1321 ). however, with a curvature source that 
captures the deviations = (^)^ — dk^.jc^^ from a constant-curvature model: 

P^'^-^P^^'ie^'ie^ ' ^'^ff-'ie^ + is^- ^^^^ 



2.3.3 'Morphed' Friedmann cosmologies 



In the above-introduced framework we distinguish the averaged matter source 
on the one hand, and averaged sources due to geometrical inhomogeneities stem- 
ming from extrinsic and intrinsic curvature (kinematical backreaction terms) on 
the other. As shown above, the averaged equations can be written as standard 
Friedmann equations that are sourced by both. Thus, we have the choice to con- 
sider the averaged model as a (scale-dependent) 'standard model' with matter 
source evolving in a mean field of backreaction terms. This form of the equations 
is closest to the standard model of cosmology. It is a 'morphed' Friedmann cos- 
mology, sourced by matter and 'morphed' by a (minimally coupled) scalar field, 
the morphon field 11481 . We write (recall that we have no matter pressure source 
here): 

pSi=--{p)9+pi ; pfn=--p% , (36) 

with 

p| = ei<^| + C/g ; pi^e^4>l-U^, (37) 

where e = +1 for a standard scalar field (with positive kinetic energy), and e = 
— 1 for a phantom scalar field (with negative kinetic energy 0- Thus, in view of 

We have chosen the letter U for the potential to avoid confusion with the volume functional; 
if e is negative, a 'ghost' can formally arise on the level of an effective scalar field, although the 
underlying theory does not contain one. 
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Eq. (|32] |. we obtain the following correspondence: 



1 



Q9 = e4>l-Us, ; 



1 



(38) 



SttG 



SttG 



Inserting ( |38] | into the integrability condition ( |29l ) then implies that for ^ 



The above correspondence allows us to interpret the kinematical backreaction ef- 
fects in terms of properties of scalar field cosmologies, notably quintessence or 
phantom-quintessence scenarii that are here routed back to models of inhomo- 
geneities. Dark Energy emerges as unbalanced kinetic and potential energies due 
to structural inhomogeneitie^^ For a full-scale discussion of this correspondence 
see ESI. 

2.3.4 A note on closure assumptions 

This system of the averaged equations in the various forms introduced above does 
not close unless we specify a model for the inhomogeneities. Note that, if the 
system would close, this would mean that we solved the scalar parts of the GR 
equations in general by reducing them to a set of ordinary differential equations on 
arbitrary scales. Closure assumptions have been studied by prescribing a cosmic 
equation of state of the form pf^^ = I5(p^,a^) |35 |, [36], or by prescribing the 
backreaction terms through scaling solutions, e.g. Qi^ 0= a"^, parametrized by a 
scaling index n |48|. We shall come back to the important question of how to 
close the averaged equations later in Subsect. 14.21 



2.4 Derived dimensionless quantities 

For any quantitative discussion it is important to provide a set of dimensionless 
characteristics that arise from the above framework. 

2.4.1 The cosmic quartet 

We start by dividing the volume-averaged Hamiltonian constraint ( |23] | by the 
squared volume Hubble functional := d^/a^ introduced before. Then, ex- 

Note that the potential is not restricted to depend only on <Prj» exphcitly. An exphcit de- 
pendence on the averaged density and on other variables of the system (that can, however, be 
expressed in terms of these two variables) is generic. 

More precisely, kinematical backreaction appears as excess of kinetic energy density over 
the 'virial Isalance', cf. Eq. dsTt . while the averaged scalar curvature of space sections is directly 
proportional to the potential energy density; e.g. a void (a 'classical vacuum') with on average 
negative scalar curvature (a positive potential) can be attributed to a negative potential energy of 
a morphon field ('classical vacuum energy'). 




(39) 
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pressed through the following set of 'parameters' E3. 
I2,f :=^(p)^ ; i2f ; 12/:=-% ; i2|:=-% , (40) 

the averaged Hamiltonian constraint assumes the form of a cosmic quartet II33I 

i2„f + i2f + i2f + i2| = 1 , (41) 

showing that the solution space of an averaged inhomogeneous cosmology is 
three-dimensional in the present framework. In this set, the averaged scalar curva- 
ture parameter and the kinematical backreaction parameter are directly expressed 
through {R) ^ and Q^, respectively. In order to compare this pair of parameters 
with the 'constant-curvature parameter' that is the only curvature contribution in 
standard cosmology to interpret observational data, we can alternatively introduce 
the pair 



; ^Iv ■■= /'d^' e^;^74(0 , (42) 

a|,//| ^ 3aj^H^J,i dt' 
being related to the previous parameters by 



I2f + ^2|v = + =■■ ■ (43) 

After a little thought we see that both sides of this equality would mimick a Dark 
Energy component, Q.^ , in a Friedmannian model. Note, in particular, that it is 
not the additional backreaction parameter alone that can play this role, but it is 
the joint action with the (total) curvature parameter, or, looking to the left-hand- 
side, it is the cumulative effect acquired during the history of the backreaction 
parameter. A positive cosmological term would require this sum, or the effective 
history, respectively, to be positive. 

2.4.2 Volume state finders 

Like the volume scale factor and the volume Hubble rate Hi^, we may intro- 
duce 'parameters' for higher derivatives of the volume scale factor, e.g. the volume 
deceleration 

^^:=-?T^ = kf + 2i2|-i2f . (44) 

Cl-Si tig I 

Following ||169| [T| (see also and references therein) we may also define the 
following volume state finders involving the third derivative of the volume scale 
factor: 

^^^^ 0^ 1 ^ 2 ^ ^^^^ 

a<s> Wg 



and 



y'3 ^ 



We shall, henceforth, call these characteristics 'parameters', but the reader should keep in 
mind that these are functional on Moreover, they are dynamically coupled. 
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The above definitions are identical to those given in lll69l fT1l. however, note the fol- 
lowing obvious and subtle differences. One of the obvious differences was already 
mentioned: while the usual state finders of a global homogeneous state in the stan- 
dard model of cosmology are the same for every scale, the volume state finders 
defined above are different for different scales. The other is the fact that the vol- 
ume state finders apply to an inhomogeneous cosmology with arbitrary 3-metric, 
while the usual state finders are restricted to a FLRW metric. Besides these there 
is a more subtle difference, namely a degeneracy in the Dark Energy density pa- 
rameter: while [169. 1] denote (with obvious adaptation) 1 — Q.^ = Q.^ we have 
from the Hamiltonian constraint d^TT) — Qq + £2^, i.e. so-called X-matter 
(Dark Energy) is composed of two physically distinct components. 



2.4.3 Cosmic equation of state and Dark Energy equation of state 

We already mentioned the possibility to characterize a solution of the averaged 
equations by a cosmic equation of state pff^ = I5{p^,a^) with := p%f/p^- 
Now, we may separately discuss (i.e. without matter source) the morphon equation 
of state that plays the role of the Dark Energy equation of state p8l. 

We can express the volume state finders through this equation of state parameter 
and its first time-derivative: 

r^ = l + ^wf(l+wf)(l^Af)-^||(l-i2,f), (48) 



and 



1 ■ & 



^" = l+>v|--— ^ , (49) 

being zero for = — 1, i.e. for the case of a (scale-dependent) cosmological 
constant. As emphasized in |1691 fT||, the above expressions have the advantage 
that one can immediately infer the case of a constant Dark Energy equation of 
state, so-called quiessence models, that here correspond to scaling solutions of 
the morphon field with a constant fraction of kinetic to potential energies P8l : 

^.=^^=^,-'Q-=24±l, (50) 

£p^t -U9V9 {R}^ w|-i ' 

where the case = (no kinematical backreaction), or wf = — 1/3 (i.e. pf + 
3p0 — 0) corresponds to the 'virial condition' 

2£^„+£p^^=0, (51) 
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obeyed by the scale-dependent Friedmannian mode0- As has been aheady re- 
marked, a non-vanishing backreaction is associated with violation of 'equilib- 
rium' . Note that a morphon field does not violate energy conditions as in the case 
of a fundamental scalar field, cf. Subsect. [3.2.1l Again it is worth emphasizing that 
the above-defined equations of state are scale-dependent. 

With the help of these dimensionless parameters an inhomogeneous, anisotropic 
and scale-dependent state can be effectively characterized. 



3 Implications and further insights: 
— qualitative views on backreaction 

Having laid down a framework to characterize inhomogeneous cosmologies and 
having understood the physical nature of backreaction effects, does not entitle us 
to draw conclusions on the quantitative importance of inhomogeneities for the 
global properties of world models. It may well be that the robustness of the stan- 
dard model also withstands this challenge. A good example is provided by New- 
tonian cosmology that is our starting point for discussing the implications of the 
present framework. 



3.1 Thoughts on Newtonian cosmology and N-body simulations 

Analytical as well as numerical models for inhomogeneities are commonly studied 
within Newtonian cosmology. Essential cornerstones of our understanding of in- 
homogeneities rest on the Euclidean notion of space and corresponding Euclidean 
spatial averages. 

3.1.1 Global properties of Newtonian models 

The present framework can also be set up for the Newtonian equations and, indeed, 
at the beginning of its development the main result on global properties of New- 
tonian models was the confirmation of the FLRW cosmology as a correct model 
describing the averaged inhomogeneous variables. Technically, this result is due 
to the fact that the averaged principal invariants, encoded in Q^, are complete 
divergences on Euclidean space sections and, therefore, have to vanish on some 
scale where we impose periodic boundary conditions on the deviation fields from 
the FLRW background. The latter is a necessary requirement to obtain unique 
solutions for Newtonian models (for details see |44|). 

This point is interesting in itself, because researchers who have set up cosmo- 
logical N-body simulations did not investigate backreaction: the vanishing of the 
averaged deviations from a FLRW background is enforced by construction. The 
same remark applies to analytical models, where a homogeneous background is 
introduced with the manifest implication of coinciding with the averaged model. 



" In the case of vanishing kinematical backreaction, the scalar field is present for our definition 
of the correspondence and it models a constant-curvature term {R)cji, = 6k<^./aj^. Alternatively, 
we could associate a morphon with the deviations Wg) from the constant-curvature model only. 
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but without an explicit proof. The outcome that a FLRW cosmology indeed de- 
scribes the average of a general Newtonian cosmology can be traced back to the 
(non-trivial) property that the second principal invariant II appearing in can 
indeed be written (like the first) as a complete divergence, cf. Eq. (154b below. Since 
this is not valid in Riemannian geometry, 'global' backreaction effects - if relevant 
- entail the need of generalizing current cosmological simulations and analytical 
models. If backreaction is substantial, then current models must be considered as 
toy-models that have improved our understanding of structure formation, but are 
inapplicable in circumstances where the dynamics of geometry is a relevant issue. 
We shall learn that (i) these circumstances are those needed to route Dark Energy 
back to inhomogeneities, and (ii) at the precision level at which currently cosmo- 
logical parameters are determined, it can already be demonstrated that backreac- 
tion might potentially be a non-negligible player in the Late Universe. 

While the last point will be touched upon in Section 4, there are a number of 
more points that improve our qualitative understanding, to which we turn now. 

3.1.2 Morphological and statistical interpretation of backreaction 

The expansion law, Eq. ( 131b , is built on the rate of change of a simple morphologi- 
cal quantity, the volume content of a domain. Although functionally it depends on 
other morphological characteristics of a domain, it does not explicitly provide in- 
formation on their evolution. An evolution equation for the backreaction term 
is missing. This fact touches on the problem of closing the hierarchy of dynamical 
evolution equations mentioned in Subsect. [2.3.41 

We shall, in this subsection, provide a morphological interpretation of that 
is possible in the Newtonian framework (the following considerations substan- 
tially rely on the Euclidean geometry of space). This will improve our understand- 
ing of what Qq! actually measures, if geometry is not considered as a dynamical 
variable. We know from previous remarks that the dynamical coupling of to 
the geometry of space sections will change this picture. 

Let us focus our attention on the boundary of the spatial domain A priori, 
the location of this boundary in a non-evolving background space enjoys some 
freedom which we may constrain by saying that the boundary coincides with a 
velocity front of the fluid (hereby restricting attention to irrotational flows). This 
way we employ the Legendrian point of view of velocity fronts that is dual to the 
Lagrangian one of fluid trajectories. Let S{x^y^z^t) = s(t) define a velocity front 
at Newtonian time f , v = VS. 

Defining the unit normal vector n on the front, n = iVi'/IV^I (the sign de- 
pends on whether the domain is expanding or collapsing), the average expansion 
rate can be written as a flux integral using Gauss' theorem, 

{®)^ = ^ I ^■^d\=^ [ vdS , (52) 

with the Euclidean volume element (f'x, and the surface element da, dS = nda. 
We obtain the intuitive result that the average expansion rate is related to another 
morphological quantity of the domain, the total area of the enclosing surface: 

(0)^ = ±^ / \^S\do ■ (53) 
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The principal scalar invariants of the velocity gradient Vij =: S^tj can be trans- 
formed into complete divergences of vector fields |69|: 

I(v,-j)-0 = V-v ; II(v,j) = a)2-a2 + l02 = iv-(v(V-v)-(vV)v) ; 

IIlKy) = ^0^ + 20{a^ + ^(0^) + OijOjkOki - OijCOiCOj 
= ^V-Qv-(v(V-v)-(v-V)v)v-(v(V-v)-(v J . (54) 



(With our assumptions ft) in the above expressions vanishes identically.) 

In obtaining these expressions, the flatness of space is used essentially. Inserting 

the velocity potential and performing the spatial average, we obtain 11381 : 

(II) = f lldh= f H \VS\^do ; (55) 

{m)^ = [ md^x = ± [ G\WS\^do , (56) 

V@ Js JdS 

where H is the local mean curvature and G the local Gaussian curvature at every 
point on the 2— surface bounding the domain. jV^I = ^ equals 1, if the instrinsic 
arc-length i of the trajectories of fluid elements is used instead of the extrinsic 
Newtonian time t. The averaged invariants comprise, together with the volume, a 
complete set of morphological characteristics known as the Minkowski Function- 
als Wa of a body: 



{s):^ I d\ = V^ ■ W,{s):=\ [ do 
Js 3 Jd^ 



] r ] r Aw 

W2{s):=- Hdo ; Wi{s) := - Gdo = —x. (57) 

3 Jds 3 Jds 3 



The Euler-characteristic X determines the topology of the domain and is assumed 
to be an integral of motion (X = if the domain remains simply-connectecQ 

Thus, we have gained a morphological interpretation of the backreaction term: 
it can be entirely expressed through three of the four Minkowski Functionals: 



Qs{s) = 6 -^--^ . (58) 



The Wa ; a = 0, 1,2,3 have been introduced into cosmology in fl341 in order 
to statistically assess morphological properties of cosmic structure. Minkowski 
Functionals proved to be useful tools to also incorporate information from higher- 
order correlations, e.g., in the distribution of galaxies, galaxy clusters, density 
fields or cosmic microwave background temperature maps ( 111011 . 11121 . 11731 . 
111741 : see the review by Kerscher |108| and references therein). Related to the 



Notice that this may provide a morphological closure condition for the hierarchy of evolu- 
tion equations. 
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morphology of individual domains is the study of building blocks of large-scale 
cosmic structure II170II . II175L 

For a ball with radius R we have for the Minkowski Functionals: 

(59) 

Inserting these expressions into the backreaction term, Eq. ( |58] l. shows that 
Qc^[s) = 0, and we have proved Newton's 'Iron Sphere Theorem', i.e. the fact 
that a spherically-symmetric configuration features the expansion properties of 
a homogeneous-isotropic modeQ Moreover, we can understand now that the 
backreaction term encodes the deviations of the domain's morphology from that 
of a ball, a fact that we shall illustrate now with the help of Steiner's formula of 
integral geometry (see also 11341 ). 

Let da'^ be the surface element on the unit sphere, then (according to the Gaus- 
sian map) da = R\R2da^ is the surface element of a 2— surface with radii of cur- 
vature R\ and /?2. Moving the surface a distance e along its normal we get for the 
surface element of the parallel velocity front: 

da' = {R, + e){R2 + e)da' = Mi±£^l±M±£!^a = {l+e2H + e'G)da , 

R1R2 

(60) 

where 

+ (6,) 



2 \Ri R2 ) ^1^2 

are the mean curvature and Gaussian curvature of the front as before. 

Integrating Eq. ( 160b over the whole front we arrive at a relation between the 
total surface area A^* of the front and A^^^ of its parallel front. The gain in volume 
may then be expressed by an integral of the resulting relation with respect to £ 
(which is known as Steiner's formula defining the Minkowski Functionals of a 
(convex) body in three spatial dimensions): 



V3^=Vs<+f de'Arg , =V^ + eArg + e^ [ Hda + ^e^f Gda . 
Jo Jd9 3 Jds) 



(62) 



An important lesson that can be learned here is that the backreaction term Qr^ ob- 
viously encodes all orders of the N-point correlation functions, since the Minkow- 
ski Funktionals have this property; it is not merely a two-point term as the form 
of Qr^ as an averaged variance would suggest. In other words, a complete mea- 
surement of fluctuations must take into account that the domain is Lagrangian and 
the shape of the domain is an essential expression of the full N-point statistics of 
the matter enclosed within (For further statistical considerations of backreac- 
tion in terms of given fluctuation spectra see II47I . II109I ). Kinematically, Steiner's 
formula shows that the volume scale factor ar^, being defined through the volume 
in Eq. ( 126b , also depends on other morphological properties of ^ in the course of 
evolution. In a comoving relativistic setting, the domain ^ is frozen into the metric 
of spatial sections, so that we also understand that an evolving geometry in general 
relativity takes the role of this shape-dependence in the Newtonian framework. 

This can be shown exphcitly by using a radially-symmetric velocity field t47l . 
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3.1.3 Backreaction views originating from Newtonian cosmology and relativistic 
perturbation theory of a FLRW background 

We may place Newtonian models, but also relativistic models that suppress the 
coupling between fluctuations, encoded in Q^t, and the geometry of space sec- 
tions, into the same category: as a rule of thumb we can say that any model that 
describes fluctuations on a Euclidean 'background space' must be rejected as a 
potential candidate for a backreaction-driven cosmology. The reason is that fluc- 
tuations in those models can be subjected to periodic boundary conditions imply- 
ing a globally (on the periodicity scale) vanishing kinematical backreaction II44I . 
The very architecture of such models is simply too restrictive to account for a 
non-vanishing (Hubble-scale) being a generic property of relativistic mod- 
els. Of course, also in those models, backreaction can be investigated (a detailed 
investigation within Newtonian cosmology may be found in |47| as well as an 
application on the abundance statistics of collapsed objects |109|), but it is then 
only a rephrasing of the known cosmic variance within the standard model of 
cosmology. Nevertheless, the potential relevance of a non-vanishing backreaction 
can also be seen in Newtonian cosmology: in f471 it was found that the magnitude 
of X2| remains small throughout the evolution, being restricted to fall off to zero 
on some scale, but the indirect influence of a non-vanishing in the interior of 
the periodic box is strongly seen in the other cosmological parameters. Thus, in- 
dependent of our statement of irrelevance of the magnitude of Qq on large scales 
in Newtonian cosmology, backreaction is clearly an important player to interpret 
cosmological parameters starting at scales of galaxy surveys, and it may here be a 
key to also understand the Dark Matter problem, cf. Subsect. R.3.4[ 

We refer to the term 'quasi-Newtonian' when we think of relativistic mod- 
els that are restricted to sit locally close to a Friedmannian state, as in standard 
gauge-invariant perturbation theory U 13] , II140I . IIT38J . their average properties 
being evaluated on Euclidean space sections II139I . Although we do not refer to 
the discussion of structure on super-Hubble scales fBSll. 0141 . ifml . ifTSOl . the 
following consideration would also apply there. The integrability condition i29i . 
in essence, spells out the generic coupling of kinematical fluctuations to the evolu- 
tion of the averaged scalar curvature. Thus, the freedom taken by a generic model 
is carried by a non-vanishing Q^t (even if small) into changes of the other cos- 
mological parameters, notably the averaged scalar curvature. If that coupling is 
absent (even if is non-zero), Eq. ( 1291 ) shows that oc ^j^j {R)r^ °= a~jj^, 
i.e. the averaged curvature evolves like a constant-curvature model, and backreac- 
tion decays more rapidly than the averaged density, (p) qi °^V^^ . In other words, 
backreaction cannot be relevant today in all models that suppress this coupling 
(we shall make this more precise in the following). Therefore, as another rule of 
thumb, we may say that any (relativistic) model that evolves curvature at or in the 
vicinity of the constant-curvature model is rejected as a potential candidate for a 
backreaction-driven cosmology II48I . 

In summary. Dark Energy cannot be routed back to inhomogeneities on large 
scales in Newtonian and quasi-Newtonian models, but a careful re-interpretation 
of cosmological parameters will have nevertheless to be envisaged. 
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3.2 Qualitative picture for backreaction-driven cosmologies 

Looking at the backreaction term Qr^, the relevant positive term that could po- 
tentially drive an accelerated expansion in accord with recent indications from 
supernovae data |8|, 1641 . 11801 (see also Leibundgut and Enqvist, this volume!^ 
is the averaged variance of the rate of expansion, cf. Eq. (127b . This term, however, 
is quadratic and the averaging operation involves a division by the square of the 
volume. How can we then expect that, in an expanding Universe, such a term can 
be of any relevance at the present time? Before we give an answer to this ques- 
tion, let us introduce a criterion for a backreaction-driven cosmology that requires 
volume acceleration, i.e. we postulate high relevance of backreaction. This can be 
done with the help of the averaged equations as has been advocated by Kolb et al. 

mmm. 

3.2.1 Acceleration and energy conditions 

Let us look at the general acceleration law (1241) . and ask when we would find 
volume acceleration on a given patch of the spatial hypersurface 111 161 . II35II36I : 

3^^A-47rG(p)^ + e^ > 0. (63) 

acg 

We find that, if there is no cosmological constant, the necessary condition Qcg > 
AnG{p)^ must be satisfied on a sufficiently large scale, at least at the present 
time. This requires that Qr^ is positive, i.e. shear fluctuations are superseded by 
expansion fluctuation^ and, what is crucial, that decays less rapidly than the 
averaged density |35|. It is not obvious that this latter condition could be met in 
view of our remarks above. We conclude that backreaction has only a chance to 
be relevant in magnitude compared with the density (e.g. as defined through the 
inequality Eq. ( 1631) today), if its decay rate substantially deviates from its 'quasi- 
Newtonian' behavior and, more precisely, its decay rate must be weaker than that 
of the averaged density (or at least comparable, depending on initial data for the 
magnitude of Early Dark Energy ISTIl . 1521). 

Another model of Dark Energy is to assume the existence of a scalar field 
source, a so-called quintessence field (others are discussed in II65I ). However, 
a usual scalar field source in a Friedmannian model, attributed e.g. to phantom 
quintessence that leads to acceleration, will violate the strong energy condition 
p + 3/? >0, i.e.: 

3-^-AnG{p + 3p) = ~AnG{pH+P<t> + 'iP'f) >0 . (64) 
a 

Note, however, that the interpretation of volume acceleration in those data relies on the 
FLRW cosmology. Backreaction could be influential and could change the interpretation of as- 
tronomical data also without featuring an accelerating phase. 

From the observational point of view this property is in accord with constraints that can be 
imposed on the averaged shear fluctuations (quantitatively discussed in 1 36 1): the universe model 
can be highly isotropic in accord with strong constraints on the shear amplitude on large scales. 
For the backreaction term it is important to independently constrain the large-scale expansion 
fluctuations that are in general not necessarily proportional to large-scale density fluctuations as 
in a linear perturbation approach at a FLRW background. Note also that the time-evolution of 
an isotropic average model must not (and in this case will not) coincide with the time-evolution 
of a FLRW background. 
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In Subsect. 12.3.31 we have introduced a mean field description of kinematical 
backreaction in terms of a morphon field. For such an effective scalar field the 
strong energy condition is not violated for the true content of the Universe, that is 
ordinary dust matter. In this line it is interesting that we can identify 'violation' of 
an effective 'strong energy condition' with the acceleration condition above (c/ 
Eqs. Oall, OS): 

3 ^ = -47rG(pi| + 3p2f) = ~AnG ({p)& + pi + 3p% ) = -AnG{p)@ + , 

(65) 

which has to be positive, if the acceleration condition (|63] l is met. 
3.2.2 Curvature-fluctuation coupling 

It is clear by now that a backreaction-driven cosmology II161I must make effi- 
cient use of the genuinely relativistic effect that couples averaged extrinsic and 
intrinsic curvature invariants, as is furnished by the integrability condition ( 129b 
(or the Klein-Gordon equation ( |39t in the mean field description). While models 
that suppress the scalar field degrees of freedom attributed to backreaction (or the 
morphon field in the mean field description), and so cannot lead to an explanation 
of Dark Energy on the Hubble scale, general relativity offers a wider range of pos- 
sible cosmologies, since it is not constrained by the assumption of Euclidean or 
constant-curvature geometry and small deviations thereof. Here, it is essentially 
the requirement that the evolution of the background geometry is suppressed (nat- 
urally in Newtonian models and through 'gauge-fixing' in gauge-invariant pertur- 
bation theory), while generically the geometry is a dynamical variable and does 
not evolve independently of the perturbations. But, how can a cosmological model 
be driven away from a 'near-Friedmannian' state, if we do not already start with 
initial data away from a perturbed Friedmannian model? How does the mecha- 
nism of the coupling between geometry and matter fluctuations work, and can this 
mechanism be sufficiently effective? 



3.2.3 The 'Newtonian anchor' 

Let us guide our thoughts by the following intuitive picture. Integral properties 
of Newtonian and quasi-Newtonian models remain unchanged irrespective of 
whether fluctuations are absent or 'turned on'. Imagine a ship in a silent water 
and wind environment (homogeneous equilibrium state). Newtonian and quasi- 
Newtonian models do not allow, by construction 1991 , that the ship would move 
away as soon as water and wind become more violent. This 'Newtonian anchor' 
is lifted into the ship as soon as we allow for the coupling of fluctuations to the 
geometry of spatial hypersurfaces in the form of the averaged scalar curvature. 
It is this coupling that can potentially drive the ship away, i.e. change the inte- 
gral properties of the cosmology. Before we are going to exemplify this coupling 
mechanism, e.g. by discussing exact solutions, let us add some understanding to 
the role played by the averaged scalar curvature. 
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3.2.4 The role of curvature 

Looking at the integral of the curvature-fluctuation-coupling, Eq. (130b . we un- 
derstand that the constant-curvature of the standard model is specified by the in- 
tegration constant ke^.. This term does not play a crucial dynamical role as soon 
as backreaction is at work. Envisaging a cosmology that is driven by backreac- 
tion, we may as well dismiss this constant altogether. In such a case, the averaged 
curvature is dynamically ruled by the backreaction term and its history. Given this 
remark we must expect that the averaged scalar curvature may experience changes 
in the course of evolution (in terms of deviations from constant-curvature), as 
soon as the structure formation process injects backreaction. This picture is actu- 
ally what one needs in order to solve the coincidence problem, i.e. the observation 
that the onset of acceleration of the Universe seems to coincide with the epoch of 
structure formation. 

This mechanism can be qualitatively understood by studying scaling solutions, 
cf. Subsect. [3.3.5l which impose a direct coupling, °= (R)^. (These scaling so- 
lutions correspond to quiessence fields, Eq. ( 1501) . and have been thoroughly studied 
by many people working on quintessence (see II171II . ||65]| and references therein.) 
In the language of a morphon field, the mechanism perturbs the 'virial equilib- 
rium', Eq. ( TTSl) , such that the potential energy stored in the averaged curvature is 
released and injected into an excess of kinetic energy (kinematical backreaction). 
Thus, in this picture, positive backreaction, capable of mimicking Dark Energy, is 
fed by the global 'curvature energy reservoir' . It is clear that such a mechanism 
relies on an evolution of curvature that differs from the evolution of the constant- 
curvature part of the standard cosmology. Indeed, as we shall exemplify below, 
already a deviation term of the form Wt^ = {R)r^ — 6k^.a~^ <^ a~^ is sufficient to 
change the decay rate of from oc to oc . 

If we start with 'near-Friedmannian initial data', and no cosmological con- 
stant, then the averaged curvature must be negative today and - if we require 
the model to fully account for A - of the order of the value that we would find 
for a void-dominated Universe p481. Thus, the determination of curvature evo- 
lution, even only asymptotically |165|, [166], is key to understand backreaction. 
The difference to the concordance model is essentially that the averaged curvature 
changes from an almost negligible value at the CMB epoch to a cosmologically 
relevant negative curvature today. This is one of the direct hints to put backreaction 
onto the stage of observational cosmology, cf. Subsect. [4.3.1| 

Let us add three remarks. First, it is not at all evident that a flat Universe is 
necessarily favoured by the data throughout the evolution II98I . This latter anal- 
ysis has been performed within the framework of the standard model, and it is 
clear that in the wider framework discussed here, the problem of interpreting as- 
tronomical data is more involved. Second, it is often said that spatial curvature can 
only be relevant near Black Holes and can therefore not be substantial. Here, one 
mistakenly implies an astrophysical Black Hole, while the Schwarzschild radius 
corresponding to the matter content in a Hubble volume is of the order of the Hub- 
ble scale. As the averaged Hamiltonian constraint ( 1231) shows, the averaged scalar 
curvature is a quantitatively competitive player that could only be 'compensated' 
(and only on a specified scale) by introducing a cosmological constant. In essence, 
a cosmologically relevant curvature contribution is tiny, but this property is shared 
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by all cosmological sources. Third, even standard perturbation theory predicts a 
scaling-law for the averaged scalar curvature that substantially differs from the 
evolution of a constant-curvature model, see Subsect. [4.2[ 
(The above qualitative picture is illustrated in detail in Rasanen's review II161II .) 



3.3 Exact solutions for kinematical backreaction 

The following families of exact solutions of the averaged equations are used to 
illustrate the mechanism of a backreaction-driven cosmology. Other implications 
of these examples are discussed in [36 1. 

3.3.1 A word on the cosmological principle 

We may separate the following classes of solutions into those solutions that re- 
spect the cosmological principle and those that do not. It is therefore worth re- 
calling the assumptions behind the cosmological principle. In the literature one 
often finds a 'strong' version that demands local isotropy of the universe model. 
More realistically, however, we should define a 'weak' version that refers to the 
existence of a scale of homogeneity: we assume that there exists a scale beyond 
which all observables do no longer depend on scale. It is beyond this scale where 
the standard model is supposed to describe the Universe on average; it is simply 
unreasonable to apply this model, even on average, to smaller scales, since the 
standard, spatially flat FLRW model has an in-built scale-independence. On the 
same grounds, isotropy can only be expected on the homogeneity scale and not 
below. Accepting the existence of this scale has strong implications, one of them 
being that cosmological parameters on that scale are representative for the whole 
Universe. If this were not so, and generically we may think of, e.g. a decay of 
average characteristics with scale all the way to the diameter of the Universe as 
in a generic fractal (or multi-fractal) distribution IIIOIII . then the cosmological pa- 
rameters of the standard model would make no sense unless the scale is explicitly 
indicated. The homogeneity scale is thought to be well below the scale of the ob- 
servable Universe and within our past-lightcone. Therefore, with this assumption, 
averaging over non-causally connected regions delivers the same values as those 
already accumulated up to the homogeneity scale 11611 . 1421 . 

We are now briefly describing some exact solutions, and we mainly have in 
mind to learn about the coupling between curvature and fluctuations. 

3.3.2 Backreaction as a constant curvature or a cosmological constant 

Kinematical backreaction terms can model a constant-curvature term as is already 
evident from the integrability condition ( 1291) . Also, a cosmological constant need 
not be included into the cosmological equations, since Qs* can play this role [33J, 
PtI . II160I . and can even provide a constant exactly, as was shown in 111 161 and 
ll36l . The exact condition can be inferred from Eq. ( |24] | and ([3T| and reads: 




(66) 
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which implies = Q^{ti) = const, as the only possible solution. Such a 'cosmo- 
logical constant' installs, however, via Eq. ( 130b , a non-vanishing averaged scalar 
curvature (even for k^. = 0): 

{R)^ = ^ - ^Q^iti) ■ (67) 
This fact has interesting consequences for 'morphed' inflationary models II120I . 



3.3.3 The Universe in an out-of-equilibrium state: a fluctuating Einstein cosmos 

Following Einstein's thought to construct a globally static model, we may require 
the effective scale-factor a^ on a simply-connected 3-manifold Z without bound- 
ary to be constant on some time-interval, hence = = and Eqs. ( 1241) and 
( |23l l may be written in the form: 

My My 

= 47rG-4-A ; (R)^ = nnG-^+3A , (68) 

with the global kinematical backreaction Q^, the globally averaged scalar 3-Ricci 
curvature {R)^, and the total restmass contained in Z. 

Let us now consider the case of a vanishing cosmological constant: A = 0. The 
averaged scalar curvature is, for a non-empty Universe, always positive, and the 
balance conditions (168b replace Einstein's balance conditions that determined the 
cosmological constant in the standard homogeneous Einstein cosmos. A globally 
static inhomogeneous cosmos without a cosmological constant is conceivable and 
characterized by the cosmic equation of state: 

{R)z = ^Ql = const. = pfff = . (69) 

Eq. i69i is a simple example of a strong coupling between curvature and fluctua- 
tions. Note that, in this cosmos, the effective Schwarzschild radius is larger than 
the radius of the Universe, 

~ I. ^ I T= ~ — 2GM2; = — ^Schwarzschild , (70) 

hence confirms the cosmological relevance of curvature on the global scale L. 
The term 'out-of-equilibrium' refers to our measure of relative information en- 
tropy, cf. Subsect. [l.3.4l in the above example volume expansion cannot compete 
with information production because the volume is static, while information is 
produced (see ll36l for more details). 

Such examples of global restrictions imposed on the averaged equations do 
not refer to a specific inhomogeneous metric, but should be thought of in the spirit 
of the virial theorem that also specifies integral properties but without a guarantee 
for the existence of inhomogeneous solutions that would satisfy this condition. (In 
061 a possible stabilization mechanism of a stationarity condition by backreac- 
tion, as opposed to the global instability of the classical Einstein cosmos, has been 
discussed.) 



30 



3.3.4 Demonstration of the backreaction mechanism: 
a globally stationary inhomogeneous cosmos 

Suppose that the Universe indeed is hovering around a non-accelerating state 
on the largest scales. A wider class of models that balances the fluctuations and 
the averaged sources can be constructed by introducing globally stationary effec- 
tive cosmologies: the vanishing of the second time-derivative of the scale-factor 
would only imply = const. =: C, i.e., aj; = + C(f — f,), where the integration 
constant as is generically non-zero, e.g. the model may emerge II74I . ifTS I from a 
globally static cosmos, as := 1, or from a 'Big-Bang', if as is set to zero. In this 
respect this cosmos does not appear very different from the standard model, since 
it evolves at an effective Hubble rate Hz °= 1/f. (There are, however, substantial 
differences in the evolution of cosmological parameters, see [36], Appendix B.) 

The averaged equations deliver a dynamical coupling relation between and 
{R)z as a special case of the integrability condition ( |29|FT : 



1 . 4C^ 
3' 



d,Qz + -dt{R)^ = — . (71) 



The cosmic equation of state of the A —free stationary cosmos and its solutions 
read l35][36l: 



pfff = -3 Peff ; Qe = (72) 



1 z . ^ Qziti) 



W, = ^-^2zMM^_ (73) 
a^ a^ 

The total kinematical backreaction QeV^ = AnGM^ is a conserved quantity in this 
case. 

The stationary state tends to the static state only in the sense that, e.g. in the 
case of an expanding cosmos, the rate of expansion slows down, but the steady 
increase of the scale factor allows for a global change of the sign of the averaged 
scalar curvature. As Eq. d73]) shows, an initially positive averaged scalar curvature 
would decrease, and eventually would become negative as a result of backreac- 
tion. This may not necessarily be regarded as a signature of a global topology 
change, as a corresponding sign change in a Friedmannian model would suggest 
(see Subsect l4TT]) . 

The above two examples of globally non-accelerating universe models evi- 
dently violate the cosmological principle, while they would imply a straightfor- 
ward explanation of Dark Energy on regional (Hubble) scales: in the latter exam- 
ple the averaged scalar curvature has acquired a piece a^^ that, astonishingly, 
had a large impact on the backreaction parameter, changing its decay rate from 
a^^ to 0= a'^' , i.e. the same decay rate as that of the averaged density. This is 
certainly enough to produce sufficient 'Dark Energy' on some regional patch due 
to the presence of strong fluctuation^ [35 1. However, solutions that respect the 

^' The constant C is determined, for the normalization oj; (?,) = 1 , by : 

6C2=6A+3fe(f,-)-(i?)£(f,). 

In | 36 | a conservative estimate, based on currently discussed numbers for the cosmological 
parameters, shows that such a cosmos provides room for at least 50 Hubble volumes. 
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cosmological principle and, at the same time, satisfy observational constraints can 
also be constructed |48 1. In this latter work, scaling solutions that we shall discuss 
now, have been exploited for such a more conservative approach. 

3.3.5 The solution space explored by scaling solutions 

In 1481 a systematic classification of scaling solutions of the averaged equations 
was given. Like the averaged dust matter density (p) ^ that evolves, for a restmass 
preserving domain as (p)^^ = (p)®, a'^ , we can look at the case where also 
the backreaction term and the averaged scalar curvature obey scaling laws, 

= Q9, a% ■ {R)s = a% , (74) 

where Q^. and Rq). denote the initial values of and {R)t^, respectively. The 
integrability condition ( |29l l then immediately provides as a first scaling solution 
( Il32l , Appendix B): 

= Q3>, ■ {R)^ = /?®, a-^ . (75) 

This is the only solution with n p.ln the case n = p, we can define a coupling 
parameter (that can be chosen differently for a chosen domain of averagin g^ 
such that Q^. oc 7?^ ; the solution reads: 

fl+3r) (n + 2) 

Qrj = r (R)^ = rR^M% ; n = -^J^ ; r = , (76) 

(with r 7^ — 1 and n 7^ —6). The mean field description of backreaction. Sub- 
sect. |2331 defines a scalar field evolving in a positive potential, if Rr^. < (and 
in a negative potential if R^. > 0), and a real scalar field, if £R(gj{r + 1 /3) < 0. 
In other words, if Rc^,. < we have a priori a phantom field for r < —1/3 and a 
standard scalar field for r > —1/3; if 7?^. > 0, we have a standard scalar field for 
r < — 1 /3 and a phantom field for r > — 1 /3. 

For the scaling solutions the explicit form of the self-interaction term of the 
scalar field can be reconstructed P8||: 



(p)^,) = ^ii±^ ('(I +.)^ j (p)^, sinh^ (^v^<^^ 



(77) 

where (p)^^. is the initial averaged restmass density of dust matter, introducing a 
natural scale into the scalar field dynamics. This potential is well-known in the 
context of phenomenological quintessence models, [1691, [124|, [TSl] and 
references therein. The scaling solutions correspond to specific scalar field models 
with a constant fraction of kinetic and potential energies of the scalar field, i.e. 
with Eq. dSOll, 

i?l„ + ^^<o. = 0. (78) 

We previously discussed the case r = ('zero backreaction') for which this con- 
dition agrees with the standard scalar virial theorem. 

We turn now to an explicit discussion of these scaling solutions summarized 
in a cosmic phase diagram in Figure [T| 

For notational ease we henceforth drop the index Si and simply write r. 
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Fig. 1 This 'cosmic phase diagram', spanned by the effective volume deceleration parame- 
ter q^, Eq. ( l44t . and the effective density parameter 12®, Eq. ( I40t . is valid for all times and 
on all scales, i.e. it can be read as a diagram for the corresponding parameters 'today' on the 
scale of the observable Universe. It represents a two-dimensional subspace { A = } of the 
full solution space that would include a cosmological constant. All the scaling solutions are 
represented by straight lines passing through the Einstein-de Sitter model in the center of the di- 
agram (1/2; 1). The vertical line corresponding to (q^;l) is not associated with a solution of the 
backreaction problem; it degenerates to the Einstein-de Sitter model (1/2; 1). This line forms 
a 'mirror': inside the cone (Case E) there are solutions with I2,f > 1 that cannot be related to 
any real-valued scalar field, but are still of physical interest in the backreaction context (mod- 
els with positive averaged scalar curvature). Models with 'Friedmannian kinematics', but with 
renormalized parameters form the line r= 1/3 (for details see |48|, Appendix A). The line r = 
are models with no backreaction on which the parameter varies (scale-dependent 'Fried- 
mannian models'). Below the line r = in the 'quintessence phase' we find effective models 
with subdominant shear fluctuations (Qc^ positive, Qq negative).The line r = — 1/3 mimics a 
'Friedmannian model' with scale-dependent cosmological constant. The line below r = —1/3 
in the 'phantom quintessence phase' represents the solution inferred from SNLS data (cf. [481), 
and the point at {q^;£2^) = (— 1.03;0) locates the late-time attractor associated with this so- 
lution. Since we have no cosmological constant here, all expanding solutions in the subplane 
q'^ <0 drive the averaged variables away from the standard model featuring a backreaction- 
driven volume acceleration of effectively isotropic cosmologies that are curvature-dominated at 
late times. 
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3.3.6 Discussion ofFigureU] 



In Figure 1 we only concentrate on the two-dimensional solution space of aver- 
aged inhomogeneous cosmologies without a cosmological constant. We further 
concentrate in this discussion only on expanding universe models; the solution 
space contains also contracting models that are equally relevant if we interpret 
this figure for smaller spatial scales; (recall that we have R^. < for r > — 1 and 
>Oforr < -1). 

A phase space analysis of the scaling solutions BSll shows that the Einstein-de 
Sitter model is a saddle point for the scaling dynamics and small inhomogeneities 
with > should make the system evolve away from it. The sign of is im- 
portant: for all the models corresponding tor>Oorr< — 1, that is the cases C,D 
and E in Figure [T] which cannot produce accelerated expansion, we have Qs* < 0. 
In other words, the kinematical backreaction is dominated by shear fluctuations, 
cf. Eq. ( |27] |. This does not necessarily mean that the universe model is regionally 
(on the scale f^) anisotropic, because in these cases kinematical fluctuations decay 
rapidly. On the other hand, cases A and B that could be responsible for an accel- 
erated expansion correspond to Qr^ > and have subdominant shear fluctuations. 
Therefore, these models can be regionally almost isotropic, although kinematical 
fluctuations have strong influence. 

Moving down the cases from Case E to Case A we first have models in which 
Qe^ decays stronger than the density; equal decay rate Qr^ oc is found on 
the line r = 1/3. This situation changes for Case C where the Friedmannian 
kinematics does no longer act as an attractor: backreaction, having a decay rate 
weaker than the density, entails an averaged curvature evolution that deviates from 
a constant-curvature Friedmannian model. Case B represents the quintessence 
phase in the scalar field correspondence, in which the averaged model accelerates, 
bounded below by the line r = — 1 /3 of a constant backreaction (exactly model- 
ing a cosmological constant on a given scale). While fitting supernovae data with 
a constant negative curvature (the line r = left to the Einstein-de Sitter model) is 
not successful, we nevertheless appreciate that such Friedmannian models would 
physically mimic the instability towards a curvature-dominated phase. Deviations 
from constant-curvature carry the averaged model into the quintessence or even 
phantom quintessence regime (Case A), in which case backreaction is growing (as 
seen within the on average negatively curved space!). In Section 4, Subsect. [4.2.2l 
we shall discuss a perturbative model that features as a leading mode a decay rate 
Qs) °^ with a deviation from constant-curvature at the same rate, (/?)^ °= . 
This (conservative) model already lies in the quintessence phase of an acceler- 
ating universe model and can be located on the line r = —1/5 in between the 
constant-curvature line and the 'cosmological constant' . Thus, in this figure and 
explicitly in Figured an explanation of Dark Energy through backreaction effects 
is expressed by the expectation that a non-perturbative model would weaken the 
leading perturbative mode further; it would certainly lie below Qq, o= . We shall 
continue this discussion in the context of perturbative solutions in Subsect. [4.2.21 
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3.3.7 Explicit inhomogeneous solutions 

If we wish to specify the evolution of averaged quantities without resorting to phe- 
nomenological assumptions on the equations of state of the various ingredients, or 
on the necessarily qualitative analysis of scaling solutions, or with specific global 
assumptions, we have to specify the inhomogeneous metric 11171 . Natural first 
candidates are the spherically-symmetric Lemaitre-Tolman-Bondi (LTB) solu- 
tions that were first employed in the context of backreaction in (5T\ and (1591. 

Considerable effort has been spent on LTB solutions and, especially recently, 
relations to integral properties of averaged cosmologies have been sought. Inter- 
estingly, 114211 also found a strong coupling between averaged scalar curvature and 
kinematical backreaction, and LTB solutions also feature an additional curvature 
piece oc on some domain There are obvious shortcomings of LTB solution 
studies that consider the class of on average vanishing scalar curvature, since in 
that class also = 11471 : also here, a non-vanishing averaged curvature is cru- 
cial to study backreaction f59|. However, there is enough motivation to quantify 
the extra effect of a positive expansion variance to fit observational data ( Ii58i and 
references below). 

The value of LTB studies or studies of other highly symmetric exact solutions 
is more to be seen in the specification of observational properties such as the lu- 
minosity distance in an inhomogeneous metric lf3] l87IIT7ll78ll25ll2Tl 4|. as well as 
Enqvist (this volume). Although interesting results were obtained, especially in 
connection with the interpretation of supernova data, care must be taken when 
determining e.g. just luminosity distances, since the free LTB functions may fit 
any data ULUJ. Generally, apart from mistakes (e.g. setting the shear to zero), 
those studies sometimes confuse integral properties of a cosmological model with 
local properties (e.g. the scale factor and a local scale factor in the given met- 
ric form). The averaged equations cannot predict luminosity distances unless one 
considers averages on the lightcone, cf. Subsect. 14.3.21 (see, however, different 
strategies proposed and pursued in |1451), f22l|, and (1331), which in turn is re- 
lated to the issue of light-propagation in an inhomogeneous Universe (see 11021 . 
lfT03l . Hon, 118311 . 121, Iil30 :l, and discussion and references in |72|). A promis- 
ing strategy to exploit the LTB solution is to consider an ensemble of spherical 
regions whose initial data are constrained by a standard Cold Dark Matter power 
spectrum, and to look at the correlated average properties of the ensem blS How- 
ever, in order to avoid matching conditions that are necessarily involved for an en- 
semble of LTB solutions, a generic collapse model in the spirit of the Newtonian 
model investigated in |1091 would facilitate such a description. 

Another possibility to construct explicit inhomogeneous metrics is, of course, 
to employ perturbative, but also non-perturbative assumptions, that will be both 
discussed in Subsection 14.21 



Rasanen (priv. comm. ) is currently looking at an ensemble of spherical regions in the spher- 
ical collapse model to describe the statistical distribution of expanding and collapsing regions, 
where the statistical properties of this ensemble are given by the peak model of structure forma- 
tion for CDM. 
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4 Future theoretical and observational strategies: 
— quantitative views on backreaction 

In this section we are going to outline several strategies towards the goal of under- 
standing the quantitative importance of backreaction effects, and to device meth- 
ods of their observational interpretation. All the topics discussed below are the 
subject of work in progress. 

4.1 Global aspects 

The question of what actually determines the averaged scalar curvature is open. 
For a two-dimensional Riemannian manifold this question is answered through 
the Gauss-Bonnet theorem: the averaged scalar curvature is determined by the 
Euler-characteristic of the manifold. Hence, it is a global topological property 
rather than a certain restriction on local properties of fluctuations that determines 
the averaged scalar curvature. If such an argument would carry over to a three- 
dimensional manifold, then any local argument for an estimate of backreaction 
would obviously be off the table. (There are related thoughts and results in string 
theory that could be very helpful here.) In ongoing work |54 42] we consider the 
consequences of Perelman's work that was mentioned in Subsect. 11.2.21 There is 
no such theorem like that of Gauss and Bonnet in three dimensions, but there are 
uniformization theorems that could provide similar conclusions. For example, for 
closed inhomogeneous universe models we can apply Poincare's conjecture (now 
proven by Perelman II154II . II155II ') that any simply-connected three-dimensional 
Riemannian manifold without boundary is a homeomorph of a 3-sphere. Ongo- 
ing work concentrates on the multi-scale analysis of the curvature distribution 
and the related distribution of kinematical backreaction on cosmological hyper- 
surfaces that feature the phenomenology we observe. All these studies underline 
the relevance of topological issues for a full understanding of backreaction in rel- 
ativistic cosmology. To keep up with the developments in Riemannian geometry 
and related mathematical fields will be key to advance cosmological research. In 
this line it should be stressed that the averaged scalar curvature is only a weak de- 
scriptor for the topology in the general 3D case, and information on the sectional 
curvatures or the full Ricci tensor is required. In observational cosmology there 
are already a number of efforts, e.g. related to the observation of the topological 
structure of the Universe derived from CMB maps (for further discussion see ||36]| 
and for topology-related issues see II190III21I . 1661 . ll9lfT0]l . 11371). 

4.2 Perturbative and non-perturbative approaches to backreaction 

There is a large body of possibilities to construct a generic inhomogeneous met- 
ric. First, there is the possibility of using standard methods of perturbation theory. 
Although the equations and 'parameters' discussed in this work can live without 
introducing a background spacetime, a concrete model for the backreaction terms 
can be obtained by employing perturbation theory (preferably of the Lagrangian 
type) and, hence, a reference background must be introduced. But, the construc- 
tion idea is (i) to only model the fluctuations by perturbation theory (the term Q^) 
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and to find the final (non-perturbative) model by employing the exact framework 
of the averaged equations. Such a model is currently investigated by paraphras- 
ing the corresponding Newtonian approximation [47|. We shall outline more in 
detail below what we expect to learn from such a model. Second, we could aim 
at finding an approximate evolution equation for by (ii) closing the hierarchy 
of ordinary differential equations that involve the evolution of shear and the elec- 
tric and magnetic parts of the projected Weyl tensor. The problem of closing such 
a hierarchy of equations is often considered in the literature and various closure 
conditions are formulated (e.g., |97|). One of them, the silent universe model [.28i , 
which assumes a vanishing magnetic part of the Weyl tensor, is found to be too 
restrictive to describe a realistic inhomogeneous Universe |i77J, II179I , II191I , so 
that we need to head for closures with non-vanishing magnetic part. In this line, 
(iii) further studies of cosmic equations of state (like, e.g., the Chaplygin state 
lf89j ') are not only a clearcut way to close the averaged equations, but also a way 
to classify different solution sectors. All these models could be subjected to (iv) 
standard dynamical system's analysis to show their stability in the phase space of 
their parameters II193II188]| . 

As already remarked above, the FLRW cosmology as an averaged model is 
found to be stable in many cases, but there is an unstable sector that just lies in the 
right corner needed to explain 'Dark Energy' . In order to analyze this instability, 
we first look at perturbation theory in Lagrangian form. The following excursion 
allows us to roughly examine the possibilities provided by perturbation theory 
and to identify the unstable mode that is of interest in the Dark Energy context, 
although we do not expect such an approach to be sufficient. We shall also begin 
to investigate non-perturbative methods below. 

4.2.1 Relativistic Lagrangian perturbation theory 

The following is a shortcut to a setup that will provide insights without entering 
a detailed perturbative analysis. The idea is to generalize the Newtonian results 
on backreaction, investigated in detail in For this purpose it is enough to 
note that in a comoving and synchronous setting the electric part of the projected 
Weyl tensor is sufficient to capture the relativistic generalization of a first-order 
Lagrangian perturbation scheme in Newtonian cosmology. This latter is furnished 
by a Lagrangian set of evolution equations for a family of trajectories, sending 
an initial (Lagrangian) position vector X' to its Eulerian position vector at time f , 
x' = f(X', f) in a Euclidean embedding space. The relativistic generalization of the 
exact spatial one-forms dx' is provided by Cartan co-frame fields 77" = r\"jdx¥^ 
deforming the local exact basis dX ' . Correspondingly, the first-order Lagrangian 
perturbation solution [301 /' = a{t)X' + E,{t)P'{X'), with a{t) solving the standard 
Friedmann equations and ^ (f) a background-dependent known function of time, 
has its analog in the relativistic deformation one-form T]" = a(f )X'' + ^ {t)V"{X') 
01O61 . 1113211 . This approximation solves the 'electric part' of the projected Einstein 
equations, written for Cartan co-frame fields, to first order. This part of Einstein's 
equations, consisting of four equations for the nine co-frame coefficients 7]"^ with 



The indices (a,b,c...) are here non-coordinate indices that just count the one-forms, as 
opposed to the coordinate indices (i,j,k...). 
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determinant 

/ := det(T]';.) = ^atcS''' rin'Xk , (79) 
can be written ggl, ||38l: 

dabiiliri-] = ; ^e«foce'''n?J?';'7l =A7-47rGp,(X') . (80) 

This system of equations is the relativistic (non-Euclidean) generalization of the 
Lagrange-Newton system ( 181b below for dust matte¥^ 

5,7 4-4]=0 ; ie£„„,e'^V|,/|';/|, =A/-47rGp,(r) . (81) 

The geometrical limit that sends the non-exact Cartan forms to the exact forms df 
(implying that the metric of space is flat) reduces the system ( |80l l to the Newto- 
nian system ( 181b . demonstrating that the comoving synchronous spacetime slicing 
considered has a clearcut Newtonian limiO 



4.2.2 A non-perturbative model for backreaction and the leading mode 

Combined with the relativistic form of Zel'dovich's model 111981 , II199I , QUI , 
straightforward generalization of the results provided in |47|| yields a backreac- 
tion term that separates into its time-evolution given by ^ (f ) and the spatial de- 
pendence on the initial displacement field given by averages over the principal 
scalar invariants of the extrinsic curvature coefficients at initial time, I/, II/, III,: 

0-9 = > (^2) 

l+^(I/>^,. + ^2(II,)^,+^3^III,)^_ 

with ri := 2 (II,-)^,^. - I (I,-)|,_. = Qrj. , and 

r2:=6(III,-)5,,-^(I,.>^,(II,.>^, ; r3:=2(I,-)5,_,(III,-)5,_-^(II,-)|_. . 

The first term in the numerator is global and corresponds to the linear damping 
factor; in an Einstein-de-Sitter universe ^ ^ oc ' . The denominator of the first 
term is a volume effect, whereas the second term in brackets features the initial 
backreaction as a leading term. 

In the early stages of structure formation with ^ (f) ^ 1 we get 

« -Qs, , (83) 
a 

This (closed) system of equations was obtained in [46] for tiie case of no background source, 
in particular A = 0, and in |29 | including backgrounds of Friedmann-Lemaitre type. The func- 
tion ^ (t) is given for backgrounds including A in 1151 . A review and alternative forms of these 
equations may be found in |69|. 

A rigorous account for this Newtonian limit, employing the full set of Einstein's equations 
that includes the 'magnetic part', will be given in |49| and |38|. In a post-Newtonian setting 
the Newtonian limit leads to the Eulerian representation of the Newtonian system, while in the 
comoving setting considered here it leads to its Lagrangian representation. 
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identical to the perturbative evolution of Q^i, functionally evaluated with the lin- 
ear approximation. In the Newtonian investigation 1147 II it was found that this latter 
solution is in very good accord with the general model corresponding to ( [82l) on 
scales larger than w 300Mpc/h, which entitles us to expect that, on large scales, a 
perturbative model for Qr^ can at best moderately improve on this solution by go- 
ing to higher orders in the perturbation scheme. Since Qrg is quadratic, this mode 
appears in a relativistic second-order perturbation solution as the leading mode 
II158I . 11161 . 112311 . although this leading term is dismissed due to its property to 
be a complete divergence in a standard perturbative settinj^. Exploiting the fact 
that on large scales we only find a small deviation of the volume scale factor 
from the Friedmannian scale factor a{t) in this scheme, we may use the exact scal- 
ing solution, cf. Subsect. [3.3.51 Q|, o= as a (conservative) prototype model for 
backreaction, arising as a first leading perturbation in the vicinity of a standard 
FLRW model. The averaged scalar curvature corresponding to this scaling solu- 
tion also evolves with the same power (/?)^ °= which again is in accord with 
the leading second-order perturbative term found in II123I . 



4.2.3 Can backreaction compete with a cosmological constant? 



Let us now look at the dimensionless characteristics ( 1401) . For the perturbative scal- 
ing modes and (/?)^ discussed in the last subsection we find = —1/5X2^, 
both are growing functions of a^^, and the relevant term that can play the role of 
Dark Energy, see Eq. ( 1431) , divided by the mass density parameter, is also growing. 
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clearly demonstrating the (global) instability of the standard model. This has to be 
compared with the corresponding fraction of a cosmological constant parameter 
with respect to the density parameter. 
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(85) 



where, with the last assumption, the index of domain-dependence is redundant. 
Looking at the respective deceleration parameters. 
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we find in both models the onset of acceleration {q^^ = 1A =0) at the time when 



Q9, 
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(87) 



Notice that in our derivation of tiie large-scale behavior of a non-perturbative Lagrangian 
model, this is not the case, in agreement with the general situation in a relativistic setting. The 
backreaction term is a complete divergence only, if the initial data have this property. This latter 
is only possible for initially Euclidean geometry. 
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Although the leading second-order perturbative mode discussed here in the form 
of a scaling solution lies in the quintessence sector, cf. Figures [T] and [51 pertur- 
bation theory is restricted to a regime close to the Friedmannian state and so, 
strictly, does not allow us to follow the scaling mode further towards a curvature- 
dominated regime. However, by extrapolating the scaling behavior of the perturba- 
tive mode into this regime, its impact is in principle competitive, even if we set out 
standard initial data for Qq).: the comparison of scaling behaviors of (i) the aver- 
aged density, being a zero-order quantity in a perturbative framework, oc (ii) 
the constant-curvature, a first-order quantity (if a flat background was perturbed), 
a~^, and (iii) the backreaction terms as second-order quantities oc feature 
decay-rates that compensate the differences in their initial conditions magnitudes, 
if the volume scale factor is assumed to evolve until ag)(z = 0) w 1000 I123II . 
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Fig. 2 The unstable sector in Figure [T] that expresses the global instability of the standard 
model is shown together with the scaling behavior of the leading perturbative mode P dis- 
cussed in this and the last subsections. Again, the volume deceleration parameter is plotted 
against the effective density parameter Q,^. This scaling mode (corresponding to the coupling 
parameter r = —1/5 for the scaling index « = — 1) is shown as a dashed line. It originates 
from the Einstein-de Sitter model in the center and ends on the curvature-dominated attractor 

= 2r/(l +r) = — 1/2. This scaling solution lies in the quintessence regime, defined by the 
mean field description of a morphon field. Recall that it lies in between the line ending atq^ =0 
(models with Friedmannian kinematics with constant negative curvature) and the line ending at 

= — 1 (a morphon modeling a cosmological constant). The indicated line NP expresses our 
expectation of a non-perturbative, non-scaling solution that would fully explain Dark Energy 
today, while starting in the vicinity of the Einstein-de Sitter model. 



Thus, the expectation is that a non-perturbative treatment, allowing for an 
evolving background, would confirm our extrapolation of the perturbative mode 
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and would even produce a further weakening of the decay rates of the backreaction 
terms, eventually coming closer to the behavior of a bare cosmological constant, 
as speculated in Figure |2l Note that such a behavior, or the more extreme case 
of a growing backreaction term corresponding to a phantom quintessence in the 
scalar field correspondence, must be understood on the grounds that we are look- 
ing at the fluctuations within a negatively curved space section. In the course of 
evolution of the averaged scalar curvature, we know that the backreaction mech- 
anism draws 'potential energy' from curvature, and converts it into an excess of 
'kinetic energy' that implies the observed weakening of the decay of fluctuations. 
It is therefore misleading to think about fluctuations as evolving on a fixed back- 
ground, i.e. in 'Newtonian terms' . In this context it is worth recalling that, if the 
employed perturbative framework is 'quasi-Newtonian', then this also implies 
that backreaction terms appear as surface terms II168I . Ill 161 . II123I , demonstrat- 
ing that we are not describing fluctuations in a curved Riemannian space section in 
which case the principal scalar invariants of extrinsic curvature fluctuations cannot 
be represented through surface terms (compare Subsect. [3.1l) . 

The fact that already a perturbative mode entails departures of the averaged 
model from the standard model (a 'global' instability) means that the architecture 
of current N-body simulations and its determining parameters of the concordance 
cosmology is challenged and it might be overrestricted for the correct descrip- 
tion of the Late Universe: a (possibly indirect) impact of a few percent would 
already have severe implications for the demand of 'high-precision' cosmology. 
This statement needs consolidation in terms of quantitative considerations, an is- 
sue that is very involved and, at present, not conclusive. We shall just add a few 
remarks below. 



4.2.4 A few words on quantitative estimates of backreaction 

Based on the above-discussed scaling behavior of backreaction that is suggested 
by perturbation theory, we may discuss typical magnitudes of backreaction that 
are expected to be reached today. Since such estimates strongly rely on an ex- 
trapolation of a perturbative mode, they are merely indicative, but they give us an 
intuition of where we stand with perturbative calculations. 

First, if we naively (i.e. without investigating a sensible re-interpretation of 
observational data within the new framework) track the perturbative scaling solu- 
tion from standard Cold Dark Matter initial data on 'some large scale' of the or- 
der of the observable Universe, then the comparison of ( |84] | with dSSl l shows that 
backreaction is expected to fall short by a large amount to fully explain Dark En- 
ergy, e.g. setting Q^^ = A we obtain with w 1000, -4 • 12^^ = 2 • \Q~^Q.f° w 
0,0015, which still lies close to the perturbative regime. The initial data taken 
assume that the intial expansion fluctuation amplitude is independent and does 
not necessarily derive from density fluctuations. Estimates in the literature range 
from values (perturbative) of 0.004 for an inhomogeneity-induced A— parameter 

|[T92J up to Qq" w -0.05 0.26 (Lyman- a absorbers in the redshift range 

z G {3.8.2} II161II . which may at best be taken as an indication of a discrepancy 
between perturbative model estimates and the way of how we interpret observa- 
tional data. 
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Second, if we look at those estimates in a scale-dependent way, i.e. taking 
into account that the influence of backreaction must be compared to A on the 
observational scales at which we postulate a Dark Energy component, then the 
answer is more sensible: taking initial data for a standard Cold Dark Matter model 
from ||47]| and translating the effect on the time-history of Qq into the relativistic 

context, we would start to explain the value of Q^" by the perturbative scaling 
mode today on scales of typically below 100 Mpc, if that region is at 2-a variance 
level in the initial conditions. For a typical such region (at 1-a) we would not 
compensate A, but would talk about a significant effect in magnitude. 

The number of pitfalls in the above considerations is, however, large. A re- 
interpretation of the other cosmological parameters in terms of their scale-depen- 
dence is mandatory, especially since the indirect influence of a non-vanishing 
backreaction on the other cosmological parameters has been found to be crucial 
and actually is expected to largely outweigh the magnitude effect in Qq (compare 
the discussion in [47]). Therefore, it might not be a good idea to judge the influ- 
ence of backreaction based on the magnitude of Qq itself. We have to investigate 
realistic models beyond perturbation theory at a fixed background, before we can 
reliably discuss quantitative estimates from models. 



4.3 Issues of interpretation of backreaction within observational cosmology 
4.3.1 A first step: a quasi-Friedmannian template metric 

The particular form of the metric for an effective approximation of the inhomoge- 
neous Universe that springs to mind has been suggested and thoroughly discussed 
by Paranjape and Singh 111481 . who consider the metric form 

^g^ = -dt^+a^-fjdx'^dx^ , (88) 

with the volume scale factor on a mass-preserving compact domain ^ that 
is specified in terms of the exact kinematical equations, and a (domain-dependent) 
effective constant curvature three-metric with coefficients yfj that, as opposed to 
01481 . may also allow for a time-parametrization of the constant-curvature ap- 
pearing in yfi . The concrete form of the 3-metric coefficients we consider reads: 

where )%(f) corresponds to the (domain-dependent) constant curvature of the 
template space at time f, and dQ^ = r^{d(j)^ + sin^ {(j))d\jf^) . 

It should be emphasized that this template metric must not be a dust solution 
of Einstein's equations I 136II . II 16211 (the effective fluid of an averaged dust model 
also features a geometrical pressure). 

The reason why we wish to allow for an explicit time-dependence of the 'cur- 
vature constant' Kqi is given by the key-insight that the constant-curvature evolu- 
tion is not identical with that of the averaged 3-Ricci curvature of an inhomoge- 
neous universe model, if the degrees of freedom in inhomogeneities (kinematical 
backreaction) are taken into account, e.g. 1.32 J . 1361 . 111611 . This effective metric 



42 



provides an alternative dynamical picture to the thoughts recently advanced by Ka- 
sai II107II . who investigated the goodness of fit to supernova data for Friedmannian 
models without cosmological constant, but different curvature parameters. Thus, 
while a single standard model without cosmological constant cannot account for 
the supernova data, two such models - if applied to low- and high-redshift data 
separately - would |107|. In Iill9il we are currently investigating this model for 
the purpose of fitting supernova data. This fit must be constrained by CMB ob- 
servations, since otherwise we could not significantly distinguish the curvature 
evolution with backreaction from the constant-curvature evolution in a narrow 
range of redshifts ED, l98l, EQ). 

This form of an effective metric can be motivated on the grounds that Ricci 
flow renormalization of the average characteristics on a bumpy geometry, cf. Suh- 
sect. ll.2.2l would produce a constant-curvature slice, but only at a given instant of 
time. In general, such a flow has singularities, if the Ricci tensor is non-positive, 
and a constant-curvature model is reached only after subsequent steps of surgery 
of the manifold. However, if we assume intrinsic curvature fluctuations (not the 
averaged curvature), i.e. terms like {{R— (^)s')^)® '•° subdominant over kine- 
matical (extrinsic curvature) fluctuations, then we may assume that the actual 
inhomogeneous metric (at one instant of time!) is already close to a constant- 
curvature metric, in which case Ricci flow smoothing may be free of singularities. 
In any case, the disclaimer of using such a simple metric for e.g. calculating lu- 
minosity distances is still that we neglect the effect of inhomogeneities on light 
propagation. This issue we address now. 



4.3.2 Averaging on the lightcone 

Here, the most important step that would considerably advance the management 
of observational data, will be to investigate the averaging formalism on the light- 
cone. Such a framework is currently being constructed llSOl . It relates not only to 
all aspects of observations in terms of distances within inhomogeneous cosmolo- 
gies, but also links directly to initial data in the form of, e.g., CMB fluctuation 
amplitudes and the integrated Sachs-Wolfe effect. Relating lightcone averages to 
cosmological model averages is also possible and is in the focus of this investiga- 
tion. For example, a closed smooth lightfront would enclose a region of space that 
is characterized by the evolution of the volume scale factor employed in this re- 
port. The consequences of a quantitative importance of an integrated backreaction 
history, described through a propagating morphon along the lightcone, are obvi- 
ous. Applying generic redshift-distance relations e.g. to galaxy surveys would put 
us in the position to better understand the actual distribution of galaxies that are 
currently mapped with the help of FLRW distances. If expansion fluctuations are 
a dominant player on large-scales, we can imagine that also the galaxy density 
maps would be affected. This attempt is non-perturbative in the sense that the 
fully nonlinear optical propagation equations are averaged; quasi-Newtonian esti- 
mates may capture (on the background-defined lightcone) localized perturbation 
magnitudes II192I , but they suffer from the same restrictions as those discussed 
in Subsect. 13.1.31 i.e. the averaged curvature of the lightcone integrated over its 
full propagation history may substantially deviate from a perturbed background- 
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defined curvature. (Compare here also the remarks on metrical properties of space- 
time at the end of the following subsection.) 

4.3.3 Direct measurement of kinematical backreaction 

If we ask whether the kinematical backreaction term is observable, the an- 
swer within a Newtonian (or quasi-Newtonian) framework is straightforward: on 
the observable domain Si, Qi^ is built from invariants of the peculiar-velocity 
gradient in a Newtonian model. Ignoring geometrical fluctuations on regional 
scales may not be unrealistic to estimate this term from high-resolution maps 
of peculiar- velocities. More precisely, we need to carefully map the gradient of 
the peculiar-velocity to build the Newtonian approximation of Qr^. We so have 
to ignore the fact that in a relativistic setting Qr^ cannot be represented through 
invariants of a gradient, which is derived from a vector field. Existing catalogues 
are, however, too small and usually, for the definition of peculiar-velocities, the 
prior of a Friedmannian model is imposed, which therefore would only return the 
cosmic variance around the assumed Friedmannian background in a likely un- 
typical patch of the Universe that is statistically affected by boundary conditions 
II197I , EH. The measurement of on small scales may also provide a negative 
value, i.e. irrelevant for a direct large-scale estimate of Dark Energy, but relevant 
for a scale-dependent evaluation of Qi?- Indications for a shear-dominated 
on scales of about lOOMpc were discussed in the Newtonian analysis BtI . Two 
papers are of particular interest here: by taking the sampling anisotropics of the 
velocity field explicitly into account, Regos and Szalay |164|, already in 1989, 
reported a large effect (40 %) of the dipole and quadrupole anisotropics on the 
estimated bulk flow of an elliptical galaxy sample; around the same time, using 
the Eulerian linear approximation, Gorski |90 | already showed that the velocity 
field is significantly correlated even on scales of 100 Mpc. The measurement of 
the shear field related to weak gravitational lensing can add further information for 
backreaction on regional scales III72I . On large scales, on the other hand, we know 
several observational data that could place constraints on the value of kinematical 
variables |71 1. 'Global' bounds on Q^, where S is of the order of the CMB scale, 
can be inferred from work of Maartens et al. III28II . III29II . 

In this context, the question whether and how close our observers have to be at 
the center of a regional 'Hubble bubble', that probes the expected negative curva- 
ture region for positive backreaction, furnishes relevant observational input II185I . 
III86I . II 1 871 III25II . 13 1 . The scale of this 'reduced curvature region' likely exceeds 
scales that have been discussed in connection with peculiar-velocity catalogues. 

Another possibility is to exploit the relation of the kinematical backreaction 
term to Minkowski Functionals, as outlined in Subsect. 13.1.21 The problem here 
is to identify the boundary of the averaging region with a surface of constant 
peculiar-velocity potential. Again we need peculiar-velocity data or, alternatively, 
a model-dependent relation between iso-density and velocity potential surfaces; 
the relativistic geometrical effects are again ignored. The boundary of the aver- 
aging region plays a crucial role, since it carries higher-order correlations of the 
velocity distribution encoding the history of structure formation, and hence the 
backreaction history that was identified as the source of the general expansion law 
( 131b . Measuring Minkowski Functionals of iso-velocity potential surfaces thus di- 
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rectly mirrors the fact that Q^? is determined through all orders of the correlation 
functions. In this line it is important to point out that, even if the fluctuations in 
number density (the first moment of the galaxy distribution) and in the two-point 
correlation function (or the power spectrum, i.e. the second moment) may not be 
significant, fluctuations may show up especially in higher moments, since those 
determine the morphology of the averaging region (the phase correlations). An 
investigation of subsets from the IRAS catalogues revealed large morphological 
fluctuations up to scales of 200 Mpc that are significant on scales of the order of 
several tens of Megaparsecs, while on the scale around 10 Mpc these fluctuations 
disappeared UllOI . 11111 . This has been confirmed by a recent analysis of SDSS 
data [95], although here deviations were not so dramatic, an issue that has to be 
(and is currently) addressed with the help of a substantially improved data set. 

A direct determination of metrical properties of spacetime rather than proper- 
ties of the matter distribution from observational data furnishes a promising pro- 
gramme that relates to all the issues outlined here II94I . [126|. This programme 
relates to the fully relativistic considerations pursued here as opposed to the prior 
of a quasi-Newtonian model that usually enters into the interpretation process. 
Here it is important to realize that, irrespective of the small magnitude of the field 
strength in a weak-field situation, its derivatives may be important. If we consider 
space to be Euclidean and the gravitational field of the mass distribution to be a 
quasi-Newtonian perturbation, then we may not correctly characterize the effect 
of intrinsic curvature that is built in a highly nonlinear way from derivatives of the 
metric tensor. There are effects due to the morphological properties of the gravita- 
tional field, e.g. the volume effect being the simplest morphological characteristic 
mentioned in Subsect. 11.2.21 As Hellaby [93] showed, a volume matching of a 
Friedmannian template model to such a distribution implies an error of 10-30 per- 
cent which may be interpreted as a volume effect in a mass-preserving smoothing 
procedure due to a factor of the order of 71^/6 with which the Euclidean volume 
and the Riemannian volume of a ball differ [41 J. Such a factor cannot be regarded 
as a perturbation of 1 . Otherwise stated: the metrical properties of space could be 
very different from Euclidean in terms of the morphology (volume, shape, con- 
nectivity) of the gravitational field, not in terms of its magnitude. 

4.3.4 A common origin of Dark Energy and Dark Matter? 

Several times we have already pointed out that the scale-dependence of observ- 
ables is key to understand the cosmological parameters in the present framework. 
Viewing observational data with this additional discrimination power of a scale- 
dependent interpretation of backreaction effects, there is furthermore a link to the 
Dark Matter problem that certainly is important to be understood in relation to 
sources, i.e. Dark Matter particles, but there is also a kinematical contribution that 
may alter existing strategies of Dark Matter search. 

Concentrating on the Dark Energy problem has led us to focussing on a posi- 
tive contribution of Qi^ on large scales. However, as already mentioned above in 
the context of peculiar-velocity catalogues, the kinematical backreaction itself 
can also be negative, and a sign-change may actually happen by going to smaller 
scales. Looking at the phenomenology of large-scale structure reveals strongly 
anisotropic patterns, so that it is not implausible that on the scales of superclusters 
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of galaxies we would fincP^ a shear-dominated Qq; < 0. Thus, again as a result of 
its scale-dependence, the kinematical backreaction parameter can potentially be 
the origin of Kinematical Dark Energy, but also of Kinematical Dark Matter Oil . 

Mapping kinematical backreaction with a 'morphon field' opens further links 
to previous studies that tried to model Dark Energy and Dark Matter by a scalar 
field (|I1441|, fE\ and references to earlier work therein). Other explanations to 
unify the description of Dark Energy and that of Dark Matter may also be put 
into perspective [84|. With this in mind, the volume deceleration functional (1441) 
can change sign too, but this crucially depends on the value of the matter density 
parameter i2,f. We infer from Eq. ( |44] | that, for a small value of Q^, a smaller 
negative value of £2q is needed to obtain volume acceleration, < 0. Since this 
problem touches on a scale-dependent understanding of cosmological parameters, 
we now propose a strategy to properly address this issue. 



4.3.5 Multi-scale analysis of backreaction 

Let us discriminate different spatial scales by a suitable partitioning of space sec- 
tions. We denote by a scale larger than the homogeneity scale, say the Hubble- 
scale, by the scale of a typical void, and by a typical scale of a matter- 
dominated region (e.g. galaxy clusters) [42J. In standard cosmology we would 
require 12,^ w 1/4 including Dark Matter. Hence, in order to find volume accel- 
eration today, cf. Eq. ( 163b . we would need —^q" > 1/16. If, however, the global 
value of the matter parameter on the scale Ly^ is smaller, then also the needed 
amount of backreaction in a Hubble-domain J^f is smaller. Now, we discuss that 
it is indeed the case that the matter density parameter drops substantially at around 
the scale Lg in a cosmological slice that is volume-dominated by voids. 

We employ the averaged Hamiltonian constraint (123 b . and assume that a do- 
main as large as Jif is formed out of a union of underdense regions and a union 
of occupied overdense regions We further consider the following picture that 
complies with what we see in the present-day Universe: we require the volume 
Hubble expansion to be subdominant in matter-dominated regions and, on the 
other hand, the averaged density to be subdominant in devoid regions. In the first 
case, an expansion or contraction would negatively contribute and so would, e.g., 
enhance a negative averaged curvature, in the second case, the presence of a low 
averaged density would positively contribute. We can therefore reasonably expect 
that the following idealization of the distributions would not substantially impair 
the overall argument: we model voids with {p)g = and matter-dominated re- 
gions with =0 (corresponding to the stable clustering hypothesis). We also 
introduce a parameter for the occupied volume fraction, := V^^/Vyf, where 
denotes the total volume of the union of occupied regions that may be 
chosen more conservatively to weaken this idealization. Thus, we would have: 

(R)^ = -6Hj-Qg + 2A ■ {R),,g = -Qj^ + \6nG{p)„a+2A , (90) 



This was actually found in the Newtonian investigation |47| that, however, suffers from the 
fact that is restricted to drop to zero on the periodicity scale of the fluctuations. 
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together witH3 Hjf = {I - ?i^)Hs and (p),ir = A #(p),^. 

Consider for the moment the case where the kinematical backreaction terms 
in the above equations are negligible and that there is no cosmological constant. 
Then, we infer that the averaged scalar curvature must be negative on domains 
and positive on domains what obviously complies with what we expect. We 
form the 'global' cosmological parameters by dividing by 'regional' cos- 
mological paramters may be introduced by dividing by //|, if we wish to relate 
sources to the regionally measured Hubble parameter. The introduction of cosmo- 
logical parameters on the scale L is pathological and useless. With our assump- 
tions the matter density parameter Q.,^^ can be traced back to the average density 
in matter-dominated regions, (p),^ = A ^(fi)j/, and thus, the global density pa- 
rameter can be reconstructed out of an observed (p) ^ on the scale L j^. Therefore, 
we find a smaller value for the density parameter on the global scale, depending 
on the value of the volume fraction of occupied regions, as a consequence of the 
compensation (through conservation of the total mass) of the missing matter in the 
regions S . 

The volume fraction is a sensible quantity since it depends on the coarsen- 
ing of the distribution. We know that even in matter-dominated regions ^ the 
matter distribution in luminous matter is very spiky leaving a lot of volume to 
empty space. Whether this argument carries over to all matter depends on how 
smoothly Dark Matter is distributed. In relativistic cosmology it is crucial that, 
unlike for the mass, there is no equipartition of curvature in Riemannian space 
sections (there is more volume available in negatively curved regions than in pos- 
itive ones); therefore, Newtonian estimates always provide a conservative upper 
limit on a realistic volume fraction. It is not implausible that a realistic value for 
X^^g could be much smaller than anticipated by Newtonian simulations that em- 
ploy a fairly large coarsening scale ( I6III : other estimates give a larger value for 
the void volume fraction, see discussion and references in II161I . ||42j). 

Finally, it should be noted that a scale-dependent analysis may be performed 
for a given slicing of spacetime, as above, but we may also expose the particular 
situation of observers, who perform measurements in matter-dominated regions, 
to a refined analysis of a scale-dependent slicing. Such a picture has been recently 
advanced by Wiltshire and coworkers [194], [195], [1221, distinguishing cosmic 
from the observer's time, and this would involve considerations of spatial renor- 
malization of average characteristics that we briefly discussed in Subsect. ri.2.2[ 



4.4 A short conclusion: opening Pandora's Jar 

Let us conclude by stressing the most important issue: quantitative relevance of 
backreaction effects. Even if all these efforts would 'only' nail down an effect 
of a few percent, rather than 75 percent, these studies would have justified their 



The fact that we expect the global Hubble parameter to be slightly smaller than the one 
measured on the scale of voids could be used, of course with more refined assumptions, to ob- 
servationally determine the volume fraction A ,^. It will be these refined assumptions together 
with a scale-dependent treatment of other relevant variables that put us in the position to seri- 
ously think about an observational determination of the void volume fraction that is certainly 
one of the key-parameters of a scale-dependent cosmology. 
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quantitative importance for observational cosmology, and what is to be expected, 
would substantially improve our understanding of the Universe. 

Especially the recent efforts, spent on the backreaction problem by a fairly 
large number of researchers, added substantial qualitative understanding to the 
numerous previous efforts that were undertaken since George Ellis initiated this 
discussion in 1984 1701 (see references in 172 1 ). The issue remains unresolved 
to date: an explanation of Dark Energy along these lines is attractive, not only 
because it naturally explains the coincidence problem. From what has been said, it 
is also physically plausible, but a reliable and unambiguous estimate of the actual 
influence of these effects is lacking. This situation may change soon and for this 
to happen it requires considerable efforts, for which some possible strategies have 
been outlined in this section. 

After those results are coming in, we may face a more challenging situation 
than anticipated by the qualitative understanding that we have. For example, while 
the explanation of Dark Energy by quintessence (or phantom quintessence) still 
allows to hide the physical consequences behind a scalar field that is open for a 
number of explanations, the mapping of a scalar field to the backreaction problem, 
as in the mean field description outlined in Subsect. 12.3.31 can no longer keep a 
phenomenological status: fluctuations exist and can be measured. There are no 
free parameters, there are initial data that can be constrained. 

Despite being premature, let us speculate that the outcome is i) a confirmation 
of the qualitative picture of a backreaction-driven cosmology, but ii) a quantita- 
tive problem to reconcile this picture with the data in the sense that there is not 
enough time for the mechanism to be sufficient. In that situation we 'lost' the stan- 
dard model for a correct description of the Late Universe, and we do not reach a 
full explanation of Dark Energy - unless - we allow for initial data that are non- 
standard. This situation would in turn ask for a comprehensive understanding of 
these required initial data, hence reconsideration of inhomogeneous inflationary 
models 180 1 and their fluctuation spectrum at the exit epoch. As further discussed 
in 061 . globally inhomogeneous initial data may arise by the very same mecha- 
nism: if backreaction plays a role due to the generic coupling of fluctuations to 
intrinsic curvature in the Late Universe, then this coupling may have been effi- 
cient also in the Early Universe. Is it conceivable that the Universe evolved out of 
a spaceform with strongly positive averaged scalar curvature that, during inflation, 
acquires 'flatness' on average, but at the end leaves an imprint in the fluctuation 
spectrum as a remnant of the kinematical conversion of curvature energy? We 
opened Pandora's Jar. 

Notwithstanding, I would consider such a situation as the beginning of a fruit- 
ful development of cosmology. As previously mentioned, the issues of scale- 
dependence of observables, the priors underlying interpretations of observations, 
the large subject of Dark Matter and, of course, the issue of Dark Energy, will 
be all interlocked and ask for a comprehensive realistic treatment beyond crude 
idealizations. 
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Appendix: Averaged ADM equations for non-vanisliing lapse function 

For completeness, we here add the general Einstein equations for a specified foli- 
ation of spacetime employing lapse and shift functions according to the Arnowitt- 
Deser-Misner, short ADM formulation Q, 111781 . and discuss the resulting system 
of spatially averaged equations for vanishing shift. 



The ADM equations recalledW\ 

Let n^ be the future directed unit normal to a three-dimensional Riemannian hy- 
persurface E. The projector into E, h^y = g^iv+n^n^ , h^^n^ = , h^yh^^ — 
h^y), induces in E the 3-metric 

hij:=g^yh^,h^j . (A.l) 

Let us write 

«^=7V(-1, 0,0,0) , n^ = ^{l,^N'), (A.2) 

with the lapse function N and the shift vector components A'^' . Note that A'^ and A'^' 
determine our choice of coordinates. 

From «^ =g^ivn^ we find goo = -{N^- NjN') ; goj = Nj; gij = hij and, using local 
coordinates x' in a f = const, hypersurface E with 3— metric gij, setting —t and 
dx^ = dt, the line element becomes: 

ds^ = -{N^- NiN' )dt^+ 2Ni dx'dt + gij dx' ® dx^ 

= -N^dt"^ +gij{dx' +N'dt)®{dx' +N-'dt) . (A.3) 

Introducing the extrinsic curvature on E by 

Kij := ^n^-yh^.h^j = -m-j , (A.4) 

we obtain ths Arnowitt-Deser-Misner, short ADM equations Q, II178I , II73I : 
Energy (Hamiltonian) constraint: 

R-K'jK'. + K^ = l6nGe + 2A , e := T^,yn^ n" ; (A.5) 

Momentum ( Codazzi) constraints: 

K'.^^f-Kiij = 87zGJj , Ji:=-T^yn^h\ ; (A.6) 

Evolution equation for the first fundamental form: 

^' Notation: a semicolon denotes covariant derivative with respect to the 4-metric with sig- 
nature (—,+,+,+) (the units are such that c = 1), a double vertical slash covariant spatial 
differentiation with respect to the 3-metric, and a single slash denotes partial differentiation 
with respect to the local coordinates; greek indices run through ... 3, and latin indices through 
1 ... 3; summation over repeated indices is understood. 
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Evolution equation for the second fundamental form: 

1 = R^,+KKK - S)A - liVll^ii, + 1 {k\n\. - K)N\^, +N^K)^^,) 

-8KG{S'j + ^d'j{e-S\)) ; Stj -.^ T^yh^-h^j . (A.8) 
For the trace parts of (A2c) and (A2d) we have: 

l^lg^2g{-K+^N\,) , g:^dct{gij) ; (A.9) 

l-l.K = R + K^-47tG{3e-S\)-3A-^Nll'^\\, + ^N%\k ■ (A.IO) 

In relativistic cosmology it is often assumed that the energy-momentum tensor 
has the form of a perfect fluid T^y = Eu^Uy + ph^y. Also, it is often required that 
the fluid is irrotational; putting the shift vector field = 0, we then model all 
inhomogeneities of the fluid by the 3-metric and the lapse function. The lapse 
function is related to the fluid acceleration in the hypersurface that reduces to the 
pressure gradient in fluid-comoving gauge (see below): 

N e+p 

Notice that with this gauge choice the unit normal coincides with the 4-velocity 
and, especially, the momentum flux density in L vanishes. The total time-derivative 
operator of a tensor field ^ along integral curves of the imit normal, d/dt^ := 
n^dy^ = u^dy^ becomes 

dT N dt 

since «^=^|v = 0- Note that, although the definition of proper time is t : = / Ndt, 
the line element cannot be written in the form of the comoving gauge by mea- 
suring 'time' through proper time dx = Ndt, since dx is not an exact form in the 
case of an inhomogeneous lapse function. The exterior derivative of the proper 
time will involve a non-vanishing shift vector according to the space-dependence 
of the lapse function. Therefore, a foliation into hypersurfaces T = const, with 
simultaneously requiring = — (9aT is not possible. 



Averaged ADM equations for vanishing shift 

For vanishing shift vector, as will be our choice for the averaged equations, the 
line element reads: 

ds^ = -N^dt^ + gijdX'®dXi , (A.13) 

where we have written the local coordinates in capital letters now, as in the main 
text, to indicate that they now label comoving fluid elements. 
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We here recall the results given in Ii34il . We shall study spatial averages in a 
hypersurface defined by the choice of the in general inhomogeneous lapse function 
N in the line-element ( IA.13I) . 

We consider perfect fluid sources, i.e. energy density £ and pressure p with en- 
ergy momentum tensor T^y = eu^Uv + ph^y Restricting attention to irrotational 
flows we can, without loss of generality, write the flow's 4-velocity in the form 

u^ = y-^-,h='-±P, (A.14) 
h p 

together with the decomposition into kinematical parts of the 4-velocity gradient, 

M/i;v = ^©/l^iv + C^v + ft'^v -M^iMv , (A.15) 

where the inhomogeneous normalization of the 4-velocity gradient h is given by 
the injection energy per fluid element and unit restmass, de = help with the rest- 
mass density p 1 100|; is the rate of expansion, a^v the shear tensor. 

The existence of a scalar 4-velocity potential ,y together with the choice 
dA.HI l implies that the conservation equations T^.^ = are satisfied, but also that 
the flow has to be irrotational and that the covariant spatial gradient of vanishes 

mmm, m-. 

(O^v = V^v'^«[«;/3] = [d^a'^/h),^^ = ; (A.16) 

with the covariant time-derivative ,5^ ■.= u^d^,S^ = h. 

We now define the averaging operation in terms of Riemannian volume inte- 
gration on the hypersurfaces orthogonal to u^, restricting attention to scalar func- 
tions f (f,Z'), 

{W{tX))rj:=^ [ W{t,X')da, , (A.18) 

with the Riemannian volume element dUg := ^d^X, g := det{gij), and the vol- 
ume of an arbitrary compact domain, V^{t) := J^Jd^X; J := y/det{gij). We in- 
troduce a dimensionless scale factor via the volume (normahzed by the volume of 
the initial domain V^j): 

- (A..9, 

This means that we are only interested in the volume dynamics of the domain; 

will be a functional of the domain's shape (dictated by the metric) and po- 
sition. We require the domains to follow the flow lines, so that the total restmass 
Mg) := pJd^X contained in a given domain is conserved. Introducing the scaled 
(t-)expansion := A'^0, the rate of change of the domain's volume within the 
spatial hypersurfaces defines the rate of volume expansion and, through (IA.19b . 
an effective volume Hubble rate: 



(A.20) 
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We shall reserve the overdot for the covariant time-derivative (defined through the 
4- velocity u^): 

d „ d Id 

— ■=11''^ = -^ , (A.21) 

and we shall abbreviate the coordinate time-derivative by a prime in all follow- 
ing equations. For an arbitrary scalar field Y{t,X') we make essential use of the 
commutation rule 

{T)'r^ - (r% = (r©)^ - {T)s{0)s, , (A.22) 

or, alternatively, (T)'^ + 3H^{T) = {T' + T&)^. A simple application is the 
proof that the total restmass in a domain is conserved: let T = p, then {p)'^ + 
3H^{p)r^ = (p' + p0)g =0 according to the local conservation law p'+p0 =0. 

We now consider the scalar parts of Einstein's equations. Their evolution is 
determined by Raychaudhuri's equation and the Hamiltonian constraint (IA.5b . The 
former can be obtained by inserting ( IA.5I ) into dA.lOl ). 

e^-^e^-2a^-4nG{£ + 3p) + £/ , (A.23) 

with the rate of shear a, .= 1 /2a' jO], and the acceleration divergence si/ : = 

{N^'^ /N) I Upon averaging these two equations, we can cast the result into a com- 
pact form (to be found under the heading Corollary 2 in II34I '): 

3^ + 4nG{£,ff + 3p,fi) = ; 
6/f| - 167rGeeff = ; 

£',ff + 3H^ (Eeff + Pes) = , (A.24) 

with the following fluctuating sources: 

lenGpes ■- l6nG{p)s< -Q<^ + \{R)s>- f^s? ; (A.25) 

e := N^e and p := N^p are the scaled energy density and pressure of matter, 
respectively. The kinematical backreaction term is given by: 

Qs> := 2{N^n)r^- ^{NQ)% ; (A.26) 

it is built from the principal scalar invariants 211 =0^ ^'j^'i ^'i = ^0 of 
the extrinsic curvature. 

The averaged scaled scalar curvature and the acceleration backreaction terms read: 
{R)^:={N^R)s ■ p^:={j^-)^ + (^^e)^ , (A.28) 
with the scaled (t-)acceleration divergence £/ N^£/ = A^^ [N^'/N) m .. 
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Some comments 

With the help of these equations more general matter models can be considered 
within the kinematically averaged framework. Notably, scalar field sources and 
radiation. As for the latter it is interesting that, due to the non-commutativity of 
averaging and time-evolution, an averaged radiation cosmos is not described by 
the familiar law in the homogeneous situation. There are source terms demonstrat- 
ing that an inhomogeneous radiation cosmos is in an out-of-equilibrium state. 
An analoguous situation occurs for the dark radiation part when averaging brane 
world cosmologies fl31 . where those source terms can be written in terms of ef- 
fective Tsallis information entropies f9&l. (Note: it is straightforward to interpret 
the averaged ADM equations for vanishing shift for the choice of a tilted slic- 
ing, i.e. where the 4— velocity is not required to coincide with the normal on the 
hypersurfaces: we have to write them for the extrinsic curvature, and not for the 
expansion tensor, which (up to the sign) agree for our choice.) 
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